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Abstract. The p X p matrix version of the r-KdV hierarchy has been recently treated as the 
reduced system arising in a Drinfeld-Sokolov type Hamiltonian symmetry reduction applied to a 
Poisson submanifold in the dual of the Lie algebra gl pr ® (C [A, A -1 ]. Here a series of extensions of 
this matrix Gelfand-Dickey system is derived by means of a generalized Drinfeld-Sokolov reduction 
defined for the Lie algebra gl pr+s (g) (C [A, A -1 ] using the natural embedding gl pr C gl pr + s for 
s any positive integer. The hierarchies obtained admit a description in terms of a p X p matrix 
pseudo-differential operator comprising an r-KdV type positive part and a non-trivial negative 
part. This system has been investigated previously in the p = 1 case as a constrained KP system. 
In this paper the previous results are considerably extended and a systematic study is presented 
on the basis of the Drinfeld-Sokolov approach that has the advantage that it leads to local Poisson 
brackets and makes clear the conformal ( W-algebra) structures related to the KdV type hierarchies. 
Discrete reductions and modified versions of the extended r-KdV hierarchies are also discussed. 
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0. Introduction 



This paper is a continuation of [1] , where it was shown how the matrix Gelfand-Dickey 
hierarchy [2,3] fits into the Drinfeld-Sokolov approach [4] (see also [5,6,7,8,9]) to generalized 
KdV hierarchies. The phase space of this hierarchy is the space of p xp matrix Lax operators 



L Vir = Pd r + md r 1 H hw r -i<9 + w r , u i eC oc (S 1 ,gl p ), 



(0.1) 



where P is a diagonal constant matrix with distinct, non-zero entries. This phase space 
has two compatible Poisson brackets: the linear and quadratic matrix Gelfand-Dickey Pois- 
son brackets. The Hamiltonians generating a commuting hierarchy of bihamiltonian flows 
are obtained from the residues of the componentwise fractional powers of the p x p di- 
agonal matrix pseudo-differential operator L p r determined by diagonalizing L p r in the 
algebra of matrix pseudo-differential operators. This system arises from a Drinfeld-Sokolov 
type Hamiltonian symmetry reduction applied to a Poisson submanifold in the dual of 
the Lie algebra gl pr <8> (L [A, A -1 ] - where gl is the central extension of the loop alge- 
bra gl = C°° (S 1 , gl pr ) — endowed with the family of compatible Poisson brackets and 
commuting Hamiltonians provided by the r-matrix (AKS) construction (see e.g. [10]). The 
corresponding reduced phase space is identified with the set of first order matrix differential 
operators C PjT of the form 



-p,r 



(0.2) 



where j p>r G C°°(S 1 , gl pr ) and A p r G gl pr <8> (C [A, A 1 ] are written as r x r matrices with 
p x p matrix entries as follows: 
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(0.3) 



where V is a p x p diagonal constant matrix for which P = (— r _1 ) r ; l pr is the pr x pr 
identity matrix. The correspondence between L Pyr in (0.1) and C Pyr in (0.2) is established 
through the relation 

C p , r (il> t 1 ,il> t 2 ,...,'tl> t r ) t = & L P ^ 1 =X^ 1 , (0.4) 

where the ipi are p-component column vectors, yielding Ui = AviA r ~ l with A := — T -1 . 

The main purpose of this paper is to derive a series of extensions of the above system 
using the natural embedding of the Lie algebra gl pr into gl pr + s for any positive integer s. 
This embedding is given by writing the general element m G gl pr + s in the block form 



m 



A 
C 



B 
D 



(0.5) 
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where A e gl pr is written as an r x r matrix with px p matrix entries, D is an s x s matrix 
and B (respectively C) is an r-component column (row) vector with p x s (s x p) matrix 
entries. In particular, the image of A p r under this embedding is A p r s e <7/ pT . +s <E> ([ [A, A -1 ], 



A, 







(0.6) 



A generalized Drinfeld-Sokolov reduction based on the Lie algebra gl pr+s <8> I [A, A -1 ] will 
be defined in such a way that the matrix Gelfand-Dickey system is recovered when setting 
all fields outside the gl pr C gl pr + s block to zero. The corresponding reduced phase space 
will turn out to be the set of first order matrix differential operators C p , r ,s of the form 



£p,r,s — lpr+s<9 ~t~ jp,r,s ~\~ A 



(0.7) 



where j p , r ,s reads 
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(0.8) 



The dynamical variables encoded in j p , r ,s reduce to those in j Pjr (0.3) upon setting the 
fields (±, t« to zero. The phase space whose points are the operators £ p , riS will be shown to 
carry two compatible local Poisson brackets, which are naturally induced by the reduction 
procedure, and an infinite family of commuting local Hamiltonians defined by the local 
monodromy invariants of £ p , r ,s- One of the Poisson brackets will be identified as a classical 
W-algebra (see e.g. [11]). This phase space can be mapped into the space of p x p matrix 
pseudo-differential operators with the aid of the usual elimination procedure: 



A^V^Vi---,^^)* = & Lp,r,a^l = A^l, 



(0.9) 



which yields 



J p,r,s 



Pd r + md^ 1 + ■■■ + u r _i<9 + u r + z+(l s d + w)~ l z-, 



(0.10) 



where ttj is related to Vi as in (0.4) and z + = V Z- = (-■ In contrast to the standard 
case, the operator L P)7 . )S attached to C Ptr ,s is now not a differential operator but contains 



a non-trivial negative part, and the mapping £ P;rjS i— > L p ^ r s is not a one-to-one mapping. 
In fact, this mapping corresponds to factoring out a residual gl s symmetry of the hierarchy 
resulting from the generalized Drinfeld-Sokokov reduction, which is generated by the current 
w G C°°(S 1 , gl s ) through the W-algebra Poisson bracket. 
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Our main results are the following: 

First, we shall prove that the set of operators Lp, r)S is a Poisson submanifold in the 
space of p x p matrix pseudo-differential operators with respect to the compatible Gelfand- 
Dickey Poisson brackets, and that the mapping C p>r ^ s \— > L p ^ iS is a Poisson mapping from 
the bihamiltonian manifold obtained as the reduced phase space in the Drinfeld-Sokolov 
reduction onto this Poisson submanifold. We shall also present the explicit form of the 
Poisson brackets on the reduced phase space. 

Second, we shall show that the mapping C PyryS h- > L p , r ,s converts the commuting Hamil- 
tonians determined by the local monodromy invariants of C Pj r j8 into the Hamiltonians gen- 
erated by the residues of the componentwise fractional powers of the diagonalized form 

L Py r,s of L p r s . 

Third, we shall derive a "modified" version of the generalized KdV hierarchy carried by 
the manifold of operators C Pj r j8 , which via the mapping £ p>r)S i— > £ p , r , s corresponds to an 
interesting factorization of L p ^ yS (given in eq. (4.21)). One of the factors in this factorization 
(the factor AK in (4.21)) arises also independently in the r = 1 case of our construction, 
when there is no Drinfeld-Sokolov reduction and we are dealing with a generalized AKNS 
hierarchy. By further factorizing this "AKNS factor" we shall obtain a second modification. 

In the p = 1, s = 1 special case, the system based on the Lax operator L p ^ jS (0.10) 
has been considered in several recent papers (see [12,13,14] and references therein) from 
the point of view of constrained KP hierarchies. Specifically, the system considered in the 
literature may be obtained from (0.10) by setting w = 0. Setting w = is consistent with 
the flows of the hierarchy, but has the inconvenient feature that the resulting Dirac brackets 
turn out to be non-local for the second Hamiltonian structure. It is interesting that one 
can in this way recover the second Hamiltonian structure postulated by Oevel and Strampp 
[13], for p = s = 1, as a non-local reduction of the local second Hamiltonian structure that 
automatically results from the Hamiltonian reduction in the Drinfeld-Sokolov approach. 
The Drinfeld-Sokolov approach to these systems that we shall present leads to a systematic 
understanding and for this reason we think it is of interest in its own right. This approach 
also has clear advantages in that it leads to local Poisson brackets and it makes clear the 
conformal, W-algebra structures (see e.g. [11]) related to these hierarchies. It is worth 
noting that in the p = 1 case the quantum mechanical versions of these W-algebras have 
been recently found to have interesting applications in conformal field theory [15]. 

The organization of the paper is apparent from the table of contents. Section 1 is a 
brief review of the version of the Drinfeld-Sokolov approach that will be used. In Section 2 
the generalized Drinfeld-Sokolov reduction relevant in our case is defined and the resulting 
reduced system is analyzed in terms of convenient gauge slices. In Section 3 the residual 
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symmetries of the reduced system are pointed out. Section 4 is devoted to describing the 
mapping of the reduced system into the space of matrix pseudo-differential operators and to 
deriving the form of the Poisson brackets in terms of the reduced space variables. This map- 
ping will be considered both in the "Drinfeld-Sokolov gauge" (0.8) and in a "block-diagonal 
gauge" which gives rise to a factorization of the Lax operator L p r)S and to a modified ver- 
sion of the generalized KdV hierarchy. Section 5 contains the pseudo-differential operator 
description of the local monodromy invariants that in the first order matrix differential 
operator setting generate the natural family of commuting Hamiltonians. A few simple 
examples are presented in Section 6 as an illustration. As a byproduct from analyzing the 
examples, we shall also derive here a second modification of the hierarchy related to a fur- 
ther factorization of the pseudo-differential Lax operator. The connection of the results of 
this paper to the group theoretic classification of generalized Drinfeld-Sokolov reductions 
is further discussed in Section 7. There are two appendices, Appendix A and Appendix B, 
containing the technical details of certain proofs. 

Those readers who are more interested in the concrete description of the reduced system 
than its derivation by means of the Hamiltonian reduction could directly turn to study 
Theorem 4.4 in Section 4 and Theorem 5.2 together with Corollary 5.3 in Section 5, which 
give the form of the reduced Poisson brackets and the commuting Hamiltonians together 
with the compatible evolution equations, respectively. 

Notational conventions. Throughout the paper, N = C°°(S 1 , N) will denote the space 
of smooth loops in N for N a Lie group, a Lie algebra, a vector space, or mat(m x n): the 
algebra of m x n matrices over the field of complex numbers (T. . All algebraic operations such 
as addition, multiplication, Lie bracket, are extended to N in the standard pointwise fashion. 
For a finite dimensional vector space Vi (or for V a finite or infinite direct sum of such vector 
spaces) if we write fi'.S 1 ^ Vi it is understood that fi G V, (or so for the components of 
/ : S 1 — > V). The symbol will stand for the element of mat(m x n) containing a single 
non-zero entry 1 at the ik position; n and l n G gl n = mat(n x n) will denote the zero and 
identity matrices, respectively. The (smooth) dual of the vector space mat(m x n) will be 
identified with mat(n x m) by means of the pairing < , >: mat(n x m) x mat(m xn)^I 
given by 

< a,v >:= / tr (av) for a G mat(n x m), v G mat(m x n). (0-11) 
■Is' 



For a (suitably smooth) complex function F on mat(m x n), the functional derivative 

5F 
8u 



G mat(n x m) at u G mat(m x n) is defined by 



d SF /*27r / fijr 

—F(u + tv)\ t =o=<^,v>= / dxtr [ ———v(x) ) Vd G mat(m x n). (0.12) 
at ou J \ou(xj 
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For F depending on several arguments we shall use partial functional derivatives defined by 
extending the above formula in a natural manner. For instance, if F depends on u G gl n for 

n = n\ +n 2 and u = ( Ul1 Ul2 | with ua G mat(nj x m) then we have 4^- = I a±1 ai2 I 

\U21 U 22 ) 3 V J) 5u V«21 « 22 J 



in terms of the partial functional derivatives 4-^- = a 



j 2 ■ 



Finally, for a Lie algebra Q, we let £{Q) := Q ® (C [A, A -1 ] denote the space of Laurent 
polynomials in the spectral parameter A with coefficients in Q. The Lie bracket is extended 
to £(G) in the standard way. The reason for which £(Q) has to be carefully distinguished 
from Q = C°° (S 1 ,Q) is that x G [0, 2n] parametrizing 5" 1 has the role of the physical 
space variable (we adopt periodic boundary condition for definiteness) , while A will appear 
essentially as a bookkeeping device. 



1. A Hamiltonian reduction approach to KdV type hierarchies 



This section is a short review of the Drinfeld-Sokolov approach to KdV type hierarchies. 
From our viewpoint, the main idea of this approach is to combine Hamiltonian symmetry 
reduction with the Adler-Kostant-Symes approach to constructing commuting Hamiltonian 
flows. We concentrate on a special case that will be used in this paper. More general 
expositions can be found in [6,7,8,9]. 

Let Q be a finite dimensional Lie algebra with an invariant scalar product "tr" . (For 
most of this paper it may be assumed that Q = gl n for some n, in which case tr really is 
the standard matrix trace.) For C_ G Q arbitrarily fixed, consider the manifold 



M:={C = d + J + \C-\J eg} 



(1.1) 



of first order differential operators for which A is a free parameter, the so-called spectral 
parameter. As is well known (see e.g. [10]), Ai can be identified with a subspace of the 
dual of the Lie algebgra Q ® (C [A, A -1 ], where Q is the central extension of the loop algebra 
Q = C°°{S 1 , Q\ and there are two compatible Poisson brackets (PBs) on Ai defined by the 
following formulae. The current algebra PB, or affine Kac- Moody algebra PB, is given by 



{f,h} 2 (J) = J^tT 



5f_ Sh 
JJ'SJ 



Sf_ (Sh\' 
Jj \5J 



(1.2) 



This PB (or its reduction) is often referred to as the second PB. The first PB reads 



{f,hh(J) = - [ trC_ 



6f Sh 



(1.3) 



Notice that the first PB { , }i is minus the Lie derivative of the second PB { , } 2 along 
translations of J in the direction of C_. The Hamiltonians of special interest here are 
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provided by evaluation of the invariants ("eigenvalues") of the monodromy matrix T(J, A) 
of £, which is given by the path ordered exponential 

T(J, X) = V exp ^— dx(J(a;) + AC_)^ , (1.4) 

where a; G [0, 27r] parametrizes the space S 1 . The corresponding Hamiltonian flows commute 
as a special case of the Adler-Kostant-Symes (r- matrix) construction (see e.g. [10]). We call 
the hierarchy of these bihamiltonian flows on M. the "AKS hierarchy". This hierarchy 
is non-local in general since the invariants of the monodromy matrix (1.4) are non-local 
functionals of J. We wish to perform a Hamiltonian reduction of the AKS hierarchy to 
obtain a local hierarchy. The locality refers both to the commuting Hamiltonians and the 
reduced PBs. 

The standard method to obtain commuting local Hamiltonians relies on a perturbative 
procedure that uses some graded, semisimple element A of the Lie algebra £(G) = G ® 
(L [A, A -1 ]. The grading in which A is supposed to be homogeneous with a non-zero grade, 
say grade k > 0, is defined by the eigenspaces of a linear operator dN,H '■ £(G) — > £(G), 

d N , H = N\-^ + &dH, (1.5) 

where A is a non-zero integer and H e Q is diagonalizable with (usually) integer eigenvalues 
in the adjoint representation. Note that adH defines a grading of G-, 

G = ® l G 1 , [Gi,Gi}cG i+ i. (1.6) 

In the cases of our interest the grade of A is small (usually 1) and hence A takes the form 

A = (C+ + AC_) with some C± E Q. (1.7) 

The requirement that A is semisimple means that it defines a direct sum decomposition 

£(G) = Ker(ad A) + Im(ad A) , Ker(ad A) n Im(ad A) = {0}, (1.8) 

where the centralizer Ker(adA) of A is a subalgebra of £(Q). Having chosen A of the 
form given by (1.7), thereby defining C_ in (1.1); and having chosen also a compatible 
grading operator dw t H, the construction then involves imposing constraints on M. so that 
the constrained manifold M. c C M. consists of operators C of the form 

C = d+(j + C+) + XC-=d+j + A with j:S 1 ^^ j G i . (1.9) 

i<k 
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That is to say, in addition to the semisimple leading term A, C contains only terms of 
strictly smaller grade than the grade of A. In concrete applications there might be further 
constraints on C G A4 C , but the above grading and semisimplicity assumptions are already 
sufficient to obtain local Hamiltonians by the subsequent procedure. The crucial step is to 
transform C in (1.9) as follows: 

(d + j + A) ^ e adS (<9 + j + A) := (d + h + A), (1.10a) 

where S = J2i<o ^» an d h = J2i<k ^ k are (formal) infinite series consisting of terms that 
take their values in homogeneous subspaces in the decomposition (1.8) according to 

S : S 1 -> (Im(adA)) <0 , h : S 1 -> (Ker(ad A)) <fc , (1.10b) 

with the subscripts referring to the grading (1.5). The above grading and semisimplicity 
assumptions ensure that (1.10) can be solved recursively, grade by grade, for both S = S(j) 
and h = h(j) and the solution is given by unique differential polynomials in the components 
of j. The Hamiltonians of interest are associated to the grade larger than —k subspace in 
the centre of the subalgebra Ker(adA) C £(G) as follows. Suppose that {Xi} is a basis of 
this linear space. The corresponding Hamiltonians are defined by 

H Xi (j):= dx (Xi,h(j(x))) , (1.11) 
Jo 

where we use the canonical scalar product 

(x ' y):= i/f tr(x(A)y(A)) (L12) 

for any X, Y G i{Q). As we shall see in examples, the local functionals Hx^j) can be 
interpreted as the coefficients in an asymptotic expansion of eigenvalues of the monodromy 
matrix of C (1.9) for A ~ oo. They will inherit the property of the monodromy invariants 
that they commute among themselves with respect to the PBs (1.2-3) on M. if the restriction 
to M. c C M. is implemented by means of an appropriate Hamiltonian symmetry reduction. 

Let G be a finite dimensional Lie group corresponding to Q and let Stab(C_) C G be 
the subgroup which stabilizes the element C_ G Q appearing in the definition (1.1) of A4. 
Let Stab(C_) be the loop group based on Stab(C_). Consider the action of Stab(C_) on 
M. given by 

(<9+ J + AC_) ^g(d + J + XC^g' 1 = g(d + J)^" 1 + AC_ , V<7 G Stab(C_). (1.13) 

The action defined by (1.13) is a symmetry of the AKS hierarchy on M. That is it leaves 
invariant the compatible PBs (1.2-3) and the monodromy invariants. If {T 1 } C Q is a basis 
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of the centralizer of C_ in Q 1 [T l , C_] = 0, then the current components J 1 = tr (TV) serve 
as the generators of this symmetry (components of the appropriate momentum map) with 
respect to the second PB (1.2). The same current components are Casimir functions with 
respect to the first PB (1.3). 

The possibility to apply Hamiltonian reduction to the AKS hierarchy on M rests on 
the action (1.13) of the symmetry group Stab(C_). The aim is to perform a symmetry 
reduction using an appropriate subgroup of this group in such a way to ensure the locality 
of the reduced AKS hierarchy. In a Hamiltonian symmetry reduction the first step is to 
introduce constraints (e.g. by restricting to the inverse image of a value of the momentum 
map). To get local Hamiltonians, the constrained manifold should consist of operators 
satisfying the conditions in (1.9). Typically, the constraints also bring a gauge freedom into 
the system, which is to be factored out to obtain the reduced system. Another requirement 
on the constraints is that the reduced PBs inherited from (1.2) and (1.3) should be given by 
local, differential polynomial formulae. This is automatically ensured if i) the reduced PBs 
are given by the original PBs of the gauge invariant differential polynomials on Ai c and ii) 
these invariants form a freely generated differential ring. In practice, property ii) is satisfied 
if the gauge orbits admit a global cross section for which the components of the gauge fixed 
current (regarded as function of j in (1.9)) define a free generating set of the gauge invariant 
differential polynomials on M. c . The existence of such gauges is a very strong condition on 
the constraints and the gauge group. Gauge slices of this type have been used by Drinfeld 
and Sokolov [4], and we refer to such gauges as DS gauges. A detailed description of the 
notion of DS gauges, including various sufficient conditions for their existence, can be found 
in [11]. 

We wish to emphasize that the presence of a non-trivial gauge freedom is an important 
characteristic of the construction of "KdV type" systems. A system of "modified KdV 
type" that can be constructed without any reference to a gauge freedom [5,6,7] does not in 
general admit a KdV type system related to it by a "generalized Miura map" . The reason 
lies in the non-trivial conditions required for the existence of a DS gauge, which is needed 
in order to obtain a KdV type system and an associated Miura map. 

We end this overview by recalling that in the original Drinfeld-Sokolov case A was 
chosen to be a grade 1 regular semisimple element from the principal Heisenberg subalgebra 
of £{Q) for Q a complex simple Lie algebra [4,16]. Requiring A to be a regular semisimple 
element means by definition that the centralizer Ker(ad A) C £(Q) appearing in the decom- 
position (1.8) is an abelian subalgebra. In the original case this centralizer is the principal 
Heisenberg subalgebra (disregarding the central extension). For Q a complex simple Lie 
algebra or gl n , the graded regular semisimple elements taken from the other Heisenberg 
subalgebras (graded, maximal abelian subalgebras) of £(Q) [17] can be classified using the 
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results of Springer [18] on the regular conjugacy classes of the Weyl group of Q (see also 
ref. [19] 1 ). One advantage of choosing A to be regular is that the centre of the centralizer 
is then maximal, giving a large set of commuting local Hamiltonians. Another important 
point is that one may associate constraints to any graded regular A in such a way that 
DS gauges are available [7,8,1,9,19]. If one implements the generalized Drinfeld-Sokolov 
reduction procedure proposed in [7] using a graded semisimple but non-regular element A, 
then the existence of a DS gauge must be separately imposed as a condition on A. This 
was noted also in the earlier work [6], where essentially the same reduction procedure was 
described. In fact, the set of A's permitted by DS gauge fixing appears rather limited [1,9]. 
These developments raise the hope to eventually be able to completely classify the KdV 
type hierarchies that may be obtained in the Drinfeld-Sokolov approach. The analysis of 
the set of examples that follows is also intended to contribute to this classification (see the 
discussion in Section 7). 

2. A generalized Drinfeld-Sokolov reduction 

We wish to apply the above formalism to Q := gl n . The graded Heisenberg subalgebras 
of £(gl n ) are classified by the partitions of n [17]. We shall choose A to be a graded semisim- 
ple element of minimal positive grade from a Heisenberg subalgebra of £(gl n ) associated to 
a partition of n into "equal blocks plus singlets" 

p times s times 

n=pr + s = r + --- + r + l + --- + l for some p > 1, r > 1, s > 1. (2.1) 

It was explained in [1] that a graded regular semisimple element only exists for the partitions 
into equal blocks n = pr or equal blocks plus one singlet n = pr + 1. When continuing 
our previous study of the equal blocks case [1] with the equal blocks plus singlet case, it 
was realized that DS gauge fixing is possible also when the partition contains an arbitrary 
number of singlets (l's). This is the motivation for the study of the more general case 
s > 1 here; and actually the analysis will be the same for any s > 1. The construction 
that will be described in this section may be also viewed as a special case of the generalized 
Drinfeld-Sokolov reduction procedure proposed in [6,7]. 

Next we introduce the necessary notation. Adapted to the partition in (2.1), an element 
m G gl n will often be presented in the following (r + 1) x (r + 1) block form: 

r r r 

m= ^ e hJ ® A i,3 + ^2 ei ' r+1 ® B i + ^2 er+1 ' i ® C * + e r+l,r+l ® D, (2.2) 
i,j = l i=l i=l 

1 In ref. [19] the graded regular semisimple elements of £(Q) were classified for Q a classical Lie 
algebra or G2 using elementary methods. The results of Springer [18], which imply this classification 
as well as its extension to any untwisted or twisted affine Kac-Moody algebra based on any simple 
Lie algebra, had been overlooked in [19]. 
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where the e^j G gl r +i are the usual elementary matrices, A i} j G gl p , D G gl s , Bi G mat(px s) 
and Ci G mat(s x p). Alternatively, we may write 



A B 



m = 



S D )' 

Introduce the r x r "DS matrix" A r G £(gl r ) given by 



(2.3) 



fO 1 
: 1 





Va 





1 

oJ 



(2.4) 



Recall that the Heisenberg subalgebra, Z P: r jS C i(gl n ), associated to the partition (2.1) is 
the linear span of the elements 



(A r ) 1 ®Y P 




°J, V/GZ, Vy p = diag(yi,...,y p ), (yt G <L , % = 1, . . . ,p), (2.5) 
together with the elements 

V/GZ, Vy fl = diag(yi,...,y a ), (y k G (C , k= l,...,s). (2.6) 



o P 

o y s A* 



Z P;r)S is a graded maximal abelian subalgebra of £(gl n ) if we choose the grading defined by 



d := d r h = f\— + adH 
dX 



(2.7a) 



(see (1.5) and (1.6)), where we take 



H := 



diag (ml p , (m — l)l p , . . . , -ml p , S ) , if r = (2m + 1) odd; 

diag (mlp, (m — l)l p , . . . , — (m — l)l p , S ) if r = 2m even. 



(2.7b) 



Incidentally, note that H equals the Cartan generator I , 



I = diag 



r — 1 _ r — 3 



l p , . . . , 



r — 1 



lp5 Os I 



(2.8a) 



of the sl2 subalgebra of gl n under which the defining representation decomposes according 
to the partition (2.1) if r is odd, and H is obtained from I by adding | to the half-integral 
eigenvalues of Iq if r is even. In the latter case an alternative grading, with respect to which 
Z Py r, s C i(gl n ) is still graded, is defined by letting the grading operator be given by 



d' := d 2r ,2i = 2 + ad/ ^ , 



(2.8b) 
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instead of by (2.7). We have adopted the grading operator (2.7) since this grading will 
permit us to treat the case of even and odd r without any difference and, in fact, this choice 
of grading does not effect the final result. 

We choose a grade 1 element A := A P)T . )S (0.6) from the Heisenberg subalgebra, 



A = 



A r <g> r o 







Q °J with r:=diag (<*!,...,<*„), (diY^idjY, d^O. (2.9) 



For generic A, all of the eigenvalues of A except for the eigenvalue are distinct. The 
eigenvalue has multiplicity s. It follows that A is a regular element of £(gl n ) if and only 
if s = 1. In addition to the linear span of the elements given in (2.5), the centralizer 
Ker(adA) C t{gl n ) of A contains the algebra £(gl s ) C £(gl n ) spanned by the elements 



P 

DX l 



V/GZ, VDegl s . (2.10) 



Hence the centre of the centralizer of A is spanned by the elements given in (2.5) together 
with the set of elements of the form 

°» ls ° A ,), V i£Z . (2.11) 

Now we consider the generalized Drinfeld-Sokolov reduction of the AKS hierarchy on 
the manifold M. defined in (1.1), with C_ and C + given by writing A as A = AC_ + C+, 
i.e., from now on 

r-l 

C_ = e r ,i <g> T and C + = ^ e i,i+i ® r - ( 2 - 12 ) 

i=i 

Using the grading of the Lie algebra Q = gl n defined by the eigenvalues of ad H , 

Q = C?_( r _i) + . . . + Q-i + Go + Gi + ■ ■ ■ + Gr-i = G<o + Go + G>o, (2-13) 
we first write the current J as 

J = J <0 + J + J>o, (2.14) 

and impose the constraint 

J>o = C+, (2.15) 
which restricts the system to the submanifold A4 C C M. given by 

M c :={C = d+j + A\jeC co (S\ G<o) }• (2.16) 

Then we factorize the constrained manifold M. c by the group M of gauge transformations 
e-f acting according to 

e f : £ ^ e f £e~ f , f e C°° (S l : G<o)- (2.17) 
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The Hamiltonian interpretation of this reduction procedure is the same as in the case of the 
standard Drinfeld-Sokolov reduction for gl n . Briefly, from the point of view of the second 
PB (1.2), it is a Marsden-Weinstein type reduction with Af being the symmetry group and 
C + a character of M (which means that the constraints defining M c C M. are first class). 
From the point of view of the first PB (1.3), the reduction amounts to fixing the values of 
Casimir functions and subsequently factoring by the group of Poisson maps M. It follows 
that the compatible PBs on M. induce compatible PBs on the space of the gauge invariant 
functions on Ai c , identified as the space of functions on the reduced space M. Te d M c /J\f. 
The reduced first PB is the Lie derivative of the reduced second PB with respect to the 
one parameter group action on the reduced space induced by the following one parameter 
group action on M. c : 

C- tC_ , W G [R . (2.18) 

This follows by combining the fact that these translations commute with the action of M 
on M c with the fact that the first PB (1.3) is the Lie derivative of the second PB (1.2) with 
respect to the respective one parameter group action on M. It is clear from (2.17) that 
the invariants (eigenvalues) of the monodromy matrix of C G M. c define gauge invariant 
functions on Ai c . In conclusion, we have a local hierarchy of bihamiltonian flows on the 
reduced phase space -M re d, which is generated by the Hamiltonians provided by the local 
monodromy invariants of C defined by the procedure in (1.10) and (1.11). 

In order to describe the reduced system in more detail, one reverts to gauge slices. We 
wish to mention two important gauges. The first is the DS gauge whose gauge section is 
the manifold Mbs C M. c given by 

-Mds := {c = d + j DS + A\ 

r 

JDS = ^2 ^ ® V r-i+l + e r,r+l ® C+ + ^r+1,1 ® C- + e r +i, P +i ® W 
i=l 

with Vi G gl p , w G gl s , ( + G mat(p x s), £_ G mat(s x p) j. (2.19) 

In explicit matrix notation j'ds := jp,r,s is given in (0.8). The space A^ds is a one-to- 
one model of M c /M with the property that, when regarded as functions on M c , the 
components of the gauge fixed current j'ds = Jds(j) provide a basis of the gauge invariant 
differential polynomials on M. c . This follows from a standard argument on DS gauge fixing 
(see e.g. [20]), which relies on the grading and the non-degeneracy condition 

Ker(adC+) n £ <0 = {0}, (2.20) 

which is satisfied in our case. The reduced AKS hierarchy on the phase space -M re d — A^ds 
is a generalization of the well-known r-KdV hierarchy, as we shall see in Sections 4 and 5. 



13 



The other important gauge is what we call the "6-gauge" (also called block-diagonal 
gauge), which is defined by the following submanifold 6 C M c , 

:= { C = d + j + A ho e C°°(S\ Go) }. (2.21) 

Correspondingly to the grading operator H in (2.7), we parametrize jo as 

j = diag (e m , a, . . . , d ) + e m+ i, r+ i <g> b + e r+ i >m+ i <g> c , (2.22a) 

where t 2 -^-] is the integral part of 2 -=^-, m = [|]. The variables on are 6>; G gl p for z 7^ 
and 

a E gl p , fcemat(pxs), c G mat(s x p), d E gl s , (2.22b) 
which we collect into the matrix 

0o:=( a c b d )^9l p+s - (2-23) 

In the G-gauge, in terms of the variables 6*^, the reduced second PB becomes just the direct 
sum of free current algebra PBs given by 

for /, h smooth functions on G. In fact, the restriction of the second PB of arbitrary gauge 
invariant functions /, h on M. c to G has the form (2.24) with / = /|e, h = h\o>. This can 
be proved in the same way as Lemma 2.1 in [1], which is of course essentially the same as 
the proof found in [4] in the scalar r-KdV case. Using DS gauge fixing we obtain a local, 
differential polynomial mapping \x : G — > Air>s yielding a generalized Miura transformation 
from the modified KdV type hierarchy on G to the KdV type hierarchy on A^ds- As is 
expected from a Miura map, the inverse of is non-local and is not single valued. In 
other words, G cannot be reached by a local gauge fixing procedure and the intersection of 
G C M. c with a gauge orbit of Af in M. c is not unique. 

Remark 2.1. The reduced current algebra PB on At c /Af is known to provide an example 
of a "classical W-algebra" [11]. A basis of the gauge invariant differential polynomials in j 
(2.16) consisting of conformal tensors is furnished by the so called lowest weight gauge. In 
order to define this gauge one considers the s/2 subalgebra S = {/_, Io, 1+} C Q satisfying 

[/o,/±] = ±/±, [/+,/_] = 2/ (2.25) 

with I + := C+, where C+ is given in (2.12). The gauge section A4\. w . C A4 C of the lowest 
weight gauge (a conformally distinguished DS type gauge) is defined by 

M.w. :={C = d + j Lw . + A I j Lw . : S 1 - Ker (ad /_) }, (2.26) 
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where Ker (ad/_) is the space of lowest weight vectors of S C Q in its representation on Q. 
An analogue of this lowest weight gauge may be associated to any s/ 2 subalgebra S C 
and it has a PB algebra induced from the current algebra which is a W-algebra [21,20]. 

Remark 2.2. In the case r = 1 it is natural to define the element A to be the diagonal 
matrix A := A Pjr= i )S given by 

A p , r =i, a := A ( £ ° ) = Adiag(di,...,d p ,0,...,0), (2.27) 



which contains s zeros and distinct, non-zero di G C for i = 1, . . . ,p. This is a grade 1 
semisimple element of £(gl n ), n = p + s, with respect to the homogeneous grading. In this 
case the Drinfeld-Sokolov reduction becomes trivial, i.e., M = M c = Mr>s = ©• For this 
reason the assumption has been made so far that r > 1. All considerations in the rest of this 
paper apply to the r = 1 case too and all of the result follow through. There are interesting 
consequences for r = 1 as well as for r > 1. 

3. Residual symmetries 

We have performed a Drinfeld-Sokolov type reduction on the system M (1.1) using 
the subgroup Af (2.17) of the symmetry group Stab(C_) defined by the stabilizer of the 
element C_ in (2.12). Here we wish to point out a residual symmetry of the reduced system 
so obtained. Consider the following transformations on M: 

£^exp^ lr ® a ^ ^ £exp ^- ^ lr ® a q))' a = dia §( a i' • • • > a p)' ( 3 - la ) 
with ctj G C 00 ^ 1 , C ) for % = 1, . . . ,p, and 

£^ exp ^V 0^£exp(-(°^ with DeC°°{S\gl a ). (3.1b) 

The group generated by these transformations is a GLi x • • • x GL\ x GL S subgroup of 
the symmetry group Stab(C_) acting on M. according to (1.13). We call it the group of 
residual symmetries and denote it by Gr. It is easily verified that these transformations 
map the constrained manifold M c C M. to itself. For grading reasons, Gr C Stab(C_) 
normalizes the gauge group N C Stab (C_), which implies that the transformations in (3.1) 
induce a corresponding action of Gr on the space of gauge orbits M c /Af. By construction, 
this induced action leaves invariant the commuting Hamiltonians as well as the compatible 
PBs of the hierarchy on M Te d = M c /Af. Notice that the gauge slices © in (2.21) and M.\. w . 
in (2.26) are mapped to themselves by the transformations in (3.1). Hence in terms of these 
gauge slices the induced action of Gr is simply given by the restriction of (3.1). The gauge 
slice M-ds in (2-19) is mapped to itself by the subgroup GL S C Gr given by (3.1b). 
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Let us write the current J G gl pr+s defining C = (d + J + AC_ ) G M. in the block form 

£)• <3 ' 2) 

where J\\ G yi pr , J22 £ fl^ s e ^ c -? similarly to the matrix m in (2.2), (2.3), With respect 
to the current algebra PB (1.2), the infinitesimal generators (the momentum map) of the 
transformations in (3.1a) and (3.1b) on M. are provided by the current components $^ and 
w given by 

$i(J) :=tr(Jn(l r (8) e«)), i = l,...,p, (3.3a) 
where en G is the usual elementary matrix, and 

w(J) := J 22 , (3.3b) 

respectively. The restrictions of these current components to Ai c (by setting J = (j + C+) 
as in(2.16)) are gauge invariant. These gauge invariant current components generate the 
induced action of the group of residual symmetries Gr on Ai r ed with respect to the reduced 
second PB. Since the Hamiltonians and the compatible PBs of the hierarchy on -M r ed are 
invariant under this group, it follows that the current components $^ and w are constants 
along the flows of the hierarchy. Of course this also follows from the fact that $j and w are 
Casimir functions with respect to the first PB (1.3). 

In conclusion, the residual symmetry (3.1) may be used to perform further reductions 
on the hierarchy obtained from the generalized Drinfeld-Sokolov reduction. 

4. The Poisson brackets on the reduced phase space 

In this section we shall find a Poisson mapping from the reduced phase space A4 c /J\f, 
endowed with the compatible PBs induced by the Drinfeld-Sokolov reduction, to the space 
of pseudo-differential operators (PDOs) with gl p valued coefficients endowed with the usual 
Gelfand-Dickey Poisson brackets [2,3,22,23]. The image includes the phase space of the 
matrix r-KdV hierarchy. The mapping will be defined by means of the elimination procedure 
similarly to the r-KdV case [4,1]. In the present case the mapping will not be one-to-one, 
but we shall be able to present the explicit form of the reduced PBs on A4 c /Af nonetheless. 

Let A be the space of pseudo-differential operators with p x p matrix coefficients: 

N 

A = {L= L s d s \L s egl p , iVGZ}. (4.1) 

s= — 00 

Multiplication of matrix pseudo-differential operators is defined in the usual way, i.e., the 
product rule is given by matrix multiplication together with the formulae 

dd' 1 = 1 and dF = Fd + F' for F G gY p , (4.2a) 
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which engender the formula 

oo 

d- 1 F = ^2(-l) i F^d- i - 1 . 

i=0 

The Adler trace [22], Tr : A -> (C , is given by 

Tr L := / tr res(L) = J tr L_i 



(4.2b) 



(4.3) 



where tr is the ordinary matrix trace. Let P± be the projectors on A onto the subalgebras 



N 



A + :={L = J2 L * dS ^ A-:={L= ]T L S «9 S }, 



(4.4) 



s=0 



respectively. Put L± := P±(L). The space .4. is a bihamiltonian manifold. For /, h smooth 
functions on A, the quadratic (second) Gelfand-Dickey PB is given by 



{f,h}^(L) = Tr 



5h L ( 5f L\ -L 6h ( L 6f 



(4.5) 



where the gradient |£ G A of / at L G .4 is defined by 



~nf(L + tA)\ t =o = Tr A-^r 



eft 



5/ 



VA g A 



(4.6) 



The Lie derivative of the quadratic bracket (4.5) with respect to the one parameter group 
of translations 

L^(L + rl p ), VtgIR, (4.7) 
is the linear (first) Gelfand-Dickey PB: 



{f,h}W(L) = Tr I L 



5L 



+ 



5h\ 
Jl) + 



-L 



5L 



Sh 
6L 



(4.8) 



which is compatible (coordinated) with the quadratic PB. 

In order to relate the Drinfeld-Sokolov reduction to the above formalism, consider the 
linear problem for £ G Ai c (2.16): 

£V = 0, V=«V&---, (4-9) 

where ip is a (pr + s) -component column vector consisting of the p-component column 
vectors ifti (i = 1, . . . , r) and the s-component column vector <fi. This system of equations 
is covariant under the gauge transformation (2.17) accompanied by the transformation 



V^e'v, feC°°(S\g <0 ). 
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(4.10) 



Observe that the component ipi is invariant under (4.10). This implies that if we derive 
from (4.9) an equation on ifti, then the operator entering that equation will be a gauge 
invariant function on M. c . The desired operator can be derived by the usual elimination 
procedure. This is particularly simple in the DS gauge (2.19), where the explicit form of 
(4.9) is 

dVi + r^2 = o, 



<9Vv-i + r?/v = o, ( 411 ) 

r 

dlp r + ^ Vitpr+l-i + (+<P = -AT^i, 
i=l 

(d + w)(f> + C-ipi = 0, 

where V is the diagonal matrix given in (2.9). If we formally solve for tpi (i = 2, . . . , r) and 
for cj) m terms of ipi using the first (r — 1) and the last equations in (4.11), respectively, and 
re-insert the result into the remaining equation, we obtain 

L^i = AV>i, (4.12) 

where 

L = A r d r + md^ 1 + u 2 d r ~ 2 + ■■■ + Ur-xd + u r + z+(l s d + w^z-, (4.13) 
with A = — r -1 and the variables being related to those in (2.19) by 

u i = Av i A r - i , z+ = -A( + , z- = C- (4.14) 



Let M C A denote the manifold of "Lax operators" L of the form (4.13). As discussed 
above, the elimination procedure gives rise to a mapping 

7T : M c -> M, 7T (£) = L, (4.15a) 

which is constant along the gauge orbits in M c . Thus we have a corresponding induced 
mapping 

7T : M c /M -> M. (4.15b) 

Observe that L in (4.13) contains only quadratic combinations of the fields z± that 
parametrize the manifold A^ds — M. c jM (2.19), but does not contain these fields in a 
linear manner. This shows that (unlike in the usual r-KdV case) the mapping 7f is not 
one-to-one. The reason for this lies in the fact that the action of the group GL S on M. c jM 
defined by (3.1b) is a symmetry of the mapping 7f , i.e., every GL S orbit in A4 C /M is mapped 
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to a single point. To understand this, observe that the original linear problem (4.9) is co- 
variant not only with respect to the gauge group M but also with respect to the group of 
residual symmetries Gr acting according to the formulae in (3.1a) and (3.1b) complemented 
with the formulae 

♦ H-pff 1 ', 9 " ))* (4.16a) 



* « exp I ( V nil*- (4.16b) 



and 

D 

Since the component ipi is invariant under the GL S action (4.16b), the operator L = tt(C) 
entering (4.12) must be also invariant with respect to the GL S action (3.1b). One may 
directly verify that this is the case as follows. If Vi, C±, w are the coordinates of a point 
C G .Mds — M c /M, then the coordinates Vi, (±, w of the GL S transformed point L G .Mds 
turn out to be 

Vi = Vi, C+=C+e" D , C_=e D C-, w = e D we~ D + e D (e- D )' . (4.17a) 

Hence 

(+(l s d + w)- 1 ^ = (+(l s d + w)- 1 ^-. (4.17b) 
Therefore the same Lax operator L is attached to C and £, 

ifOC) = *(£), (4.18) 

showing that 7f (4.15b) maps every GL S orbit in M. c jM to a single PDO as claimed above. 
Correspondingly, for the infinitesimal generator w (3.3b) of the GL S symmetry and for any 
function T of L, we have 

{w,F} 2 = 0. (4.19) 

In particular, the expansion of L in powers of d contains only such AA-invariant differential 
polynomials in the components of C G M. c which commute with w under the second PB. 

The elimination procedure may be performed on the linear problem (4.9) in the ©-gauge 
(2.21) analogously as was done above in the DS gauge (2.19). We obtain 

L e ^i = AV>i (4.20) 

with the factorized Lax operator 

L e = (A(d + fl.^])) • • • (A(d + 0_O) (AK) (A(d + 9,)) ■ ■ ■ (A(d + 9 m )) , (4.21) 
where A = — r _1 and the operator K is given by 

K={l p d + a)-b{l s d + d)- 1 c. (4.22) 
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Here 6>j, a,b,c,d are the fields parametrizing the 6-gauge according to (2.22). Since V>i is 
gauge invariant, the operators L and Lq attached to such points of Mds and 6 that lie 
on the same gauge orbit are equal: That is Lq (4.21) is a factorized form of L (4.13). In 
particular, 

M = M e , (4.23) 

where Mq := {Lq} is the set of operators Lq (4.21). Note that if r = 2m is even, there 
appear (m — 1) factors before (AK) and m factors after (AK) in (4.21). 

Remark 4-L Following Section 2, we can define an alternative version of the reduction 
procedure using the same A (2.9), but replacing H in the grading operator (2.7) by 

1 r - 1 

H k :=J + &diag (l pr ,0 s ) with any k = 0, ±-, ±1, . . . , ±— ^— , (4.24) 

where Iq is the sl 2 generator (2.8a). In fact, the reduced phase space would be the same 
for any k since Mbs would still be a global gauge slice. Different ©-type gauges would 
arise, giving rise to factorizations of the Lax operator L similar to (4.21), but with a factor 
like K in (4.22) appearing in a different position. The associated modified KdV systems 
would all be isomorphic by permutations of the variables, except for the possibility of a 
further factorization of the factor K, which will be discussed in Section 6. Modifications of 
the 7o-g ra ding like those given by the above were introduced in [20] in connection with 
gradings compatible with a given sl 2 embedding. This method was also used in [24], where 
the factorization (4.21) (but not the results below) was derived in the case p = s = 1. 

We now consider the relationship between the set Mk of operators K of the form (4.22) 
and the space Go := gl p+s - Parametrizing the general element 9q G ©o as 6>o = ^ 
like in (2.23), we have the mapping 

V :e ^M K , r}{e Q ):=l p d + a-b{l s d + d)- 1 c. (4.25) 

The natural PB on the space ©o is given by the appropriate term in (2.24), 

for /, h smooth functions on Go- We have the following result. 

Proposition 4.1. The set Mk C A of operators K of the form (4-22) is a Poisson 
submanifold of A with respect to the quadratic G elf and- Dickey PB (4-5). The mapping r] 
(4-25) is a Poisson mapping with respect to the free current algebra PB (4-26) on Go and 
the quadratic G elf and- Dickey PB on Mk- 
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Proof. Let / and h be arbitrary (smooth) functions on A. The statement of the proposition 
is equivalent to the equality 

{f,h}^or J = {fo V ,hor ] }, (4.27) 

where the l.h.s. is determined by the quadratic Gelfand-Dickey PB (4.5) and the r.h.s. 
is determined by the free current algebra PB (4.26). This equality can be verified by 
a straightforward computation. Since the computation is rather long, we relegated it to 
Appendix A. Q.E.D. 

Thanks to Proposition 4.1, we are now ready to describe the relationship between the 
PBs on M. c /M induced from the PBs (1.2-3) on M by the Drinfeld-Sokolov reduction and 
the Gelfand-Dickey PBs on M C A. 

Theorem 4.2. The set M C A of operators L (4-13) is a Poisson submanifold with respect 
to the quadratic Gelfand-Dickey PB (4-5). The mapping rr : A4 c /Af — > M (4- 15b) defined 
by the elimination procedure is a Poisson mapping, where M is endowed with the quadratic 
Gelfand-Dickey PB and M. c /N is endowed with the reduced second PB resulting from the 
current algebra PB (1.2) on M. by means of the Drinfeld-Sokolov reduction. 

Proof. We have seen that Mx C A is a Poisson submanifold with respect to the PB (4.5). 
It is well-known (see [1] Section 2.2.) that the other factors {A(<9 + 6i) \ 6i G gl p } C A 
appearing in Lq> (4.21) are also Poisson submanifolds with respect to the PB (4.5), which 
coincides with the free current algebra on these submanifolds: 

^> M w-/. b (*[s-a-s(S)> 

for /, h smooth functions of 0$. Recall the "product property" of the quadratic bracket 
according to which the product of Poisson submanifolds is also a Poisson submanifold. The 
statement of the theorem follows from this on account of (4.23) and the fact that in the 
©-gauge (2.21) the reduced second PB is given by the current algebra (2.24). Q.E.D. 

So far we have dealt with the reduced second PB on the bihamiltonian manifold M c /Af. 
Remember that the reduced first PB on M c /Af : which results from the PB (1.3) on M, 
is the Lie derivative of the reduced second PB with respect to the infinitesimal generator 
of the one parameter group action on A4 c /Af induced by the one parameter group action 
(2.18) on M c . The manifold M of Lax operators (4.13) is also a bihamiltonian manifold 
since the linear Gelfand-Dickey PB {■, (4.8) on A can be restricted to M C A. To see 
this recall that the linear PB {■, (4.8) on A is the Lie derivative of the quadratic PB 
{■, -}( 2 ) (4.5) on A with respect to the infinitesimal generator of the one parameter group 
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action (4.7) on A and notice that this group maps McAto itself. This together with the 
first statement of Theorem 4.2 implies that M C A is in fact a Poisson submanifold also 
with respect to the linear PB (4.8). Theorem 4.2 and Theorem 4.3 below state that the 
mapping tt (4.15b) is a Poisson mapping of bihamiltonian manifolds. 

Theorem 4.3. The manifold M C A is a Poisson submanifold with respect to the linear 
G elf and- Dickey PB (4-8) on A and the mapping tt : Ai c /Af — > M (4- 15b) is a Poisson 
mapping with respect to the reduced first PB (1.3) on M. c jM and the linear G elf and- Dickey 
PB (4.8) on M. 

Proof. Using the identification M. c /M ~ M.dq, it is enough to show that the mapping 
tt : .Mds M (4.15) intertwines the one parameter group action (2.18) on *Mds an d 
the one parameter group action (4.7) on M. By (4.14), this follows from the elimination 
procedure that converts £ e Mds appearing in (4.9) into L e M appearing in (4.12). 
Q.E.D. 

The above results may be used to determine the reduced first and second PBs between 
such functions on M re d = M c /N which are of the form T o tj , Hon with functions JF, H 
on M . We now wish to present the explicit formulae for the PBs of arbitrary functions on 
•Mred- We can parametrize M. re & — -Mds by the variables Ui, z±,w or equivalently by the 
variables £,z±,w where i is the positive part of L, 

r 

1 = A r d r + J2 u i dr ~^ L = £ + z + (l s d + w)- 1 z-. (4.29) 
i=i 

With the aid of the usual functional derivatives, the arbitrary variation SH of a function H 
on .Mred may be written as 




(4.30a) 



or equivalently as 

w =K£ M ) + M^ + £ fc - + *M (430b) 

using the Adler trace (4.3) and the definition 

^ = Va— ^. (4.30c) 

i=i 

We can write the reduced first and second PBs, denoted by {F, H}* (i = 1,2), of the 
arbitrary functions F, H as follows: 

{ F, H Yi = X%(F) = Tr (f K H (/)) + £ tr {^ H {z +) + jfx^-) + ^hH) , 

(4.31) 
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where is the corresponding Hamiltonian vector field associated to the function H and 
~X. l H (G) = (5G, X^) is the derivative of the function G with respect to X^. An arbitrary 
PDO A can be expanded either in the right form A = ^ fc A^d k or in the left form A = 
d k A k . Using these expansions we define 



P (A) :=A , pi(A) :=A . 
It will be also convenient to rewrite the PDO (d + w) as 

(d + w) = W~ 1 dW with W = W- 1 W, 



(4.32) 



(4.33) 



where the GL S valued function W on [R is uniquely associated to w by (4.33) and the 
condition W(0) = l s . In terms of these notations we can now write down X^. 

Theorem 4.4. The Hamiltonian vector field X^ associated to a function H on M. c /M by 
means of the reduced second PB is given by 



+ 



SH 



+ ( <^-{<) ■ w) ': 



Sz. 



- [z+id + w)- 1 ^ 



x 2 H (z + ) = p (e 



SH SH 



Sz. 



W' 1 ) W-z+ 



SH 

Sw ' 



Y 2, x SH 5H 
X H (w) = y~ z + ~ z -~ 1" 



Sz 



+ 



Sz. 



8_H_ 

Sw ' 



w 



- 



/S_H\' 

\ Sw 



The Hamiltonian vector field X^ corresponding to the reduced first PB reads 

5H~ 

+ 



X^) = 



SI 



X^(z±) 



± f , xJ r( «,) = o. 



(4.34) 



(4.35) 



Proof. On account of Theorem 4.2, in order to verify (4.34) it is enough to compute 
the Hamiltonian vector fields separately for functions on J\4 c /J\f that have the special form 
j sl tr (f±z±) or f sl tr (aw) with some matrix valued test functions f±, a. This computation 
is presented in Appendix B. After writing down the formula of the reduced second PB from 
(4.34), it is easy to compute its Lie derivative with respect to the vector field V on M. re & 
given by 

V(£) = l p , V(z±) = V(w) = 0. (4.36) 
We know from (2.18) that this gives the formula of the reduced first PB and we find 



+ / tr 



SF SH SF SH 



Sz+ Sz- Sz- Sz + 



(4.37) 
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which is equivalent to (4.35). For computational details, see Appendix B. Q.E.D 

Note that the introduction of the "integrating factor" W in the above is only a muta- 
tional trick which we used to get compact formulae. For instance, 

p (^)^ ( -x)-^(-) + e^(^), 

(4.38a) 

W -> P < «L,) = v + ±yy-^ g„) , 

where for arbitrary sxp and p x s matrix valued functions /5 and /3 on S 1 we define their 
covariant derivatives 

V{(3) := {(3' + w(S), V0) := 0' - 0w). (4.38b) 

All other terms containing W can be rewritten in terms of w in an analogous fashion. 

Let us now consider a function H that depends on £,z±,w only through the Lax 
operator L in (4.29), 

H(£,z+,z_,w) = H(L), (4.39) 
i.e., H = H o n for some function H on M. Naturally, in this case we have the equality 

5H = Tr (j^SLj . (4.40) 

Comparing (4.40) and (4.30b) using (4.29) leads to 

{jl)_ = ^ (4 ' 41a) 

up to terms that do not contribute in (4.40), and to 
SH 



^ ip »K® + )- £=*(© + *H* (44ib) 



- = -res rWz. ^ j z+W-'d-'W ] . (1.11c) 

The PDO |^ may be regarded as an arbitrary solution of these requirements. Using these 
requirements and certain identities, e.g. the identities 



L 



+ y= z+ ( d+w rw-^{wz_{^) + ) 

l ) =jP °((^0 z + w ~ i ^j w ( d + w )~ iz - 



sn 

SL 

(4.42) 

Jl 
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for H in (4.39) we can rewrite the Hamiltonian vectors fields of Theorem 4.4 as 



tf)/- l @K + 

sn 



X 2 H (z + )=P ^L—z + W- 1 jW, 

X%(z.) = -W- 1 PS(wz.^ ; L 
X.%(w) = 0, 



(4.43) 



and respectively as 



+ 



L > -FT 



sl) + 



Jl 



z+w-^j w, 

X 1 H (z.) = -W- 1 P^ 
X-hW = o. 



(4.44) 



Wz. 



6H 
~5L 



These formulae are consistent with the claims of Theorems 4.2 and 4.3. They will be used 
at the end of Section 5 to determine the evolution equations of the KdV type hierarchy that 
results from the Drinfeld-Sokolov reduction. 

Remark 4-2. The compact presentation of the formulae for the PBs in (4.43), (4.44), (4.34) 
and (4.35) was suggested by the formulae in [13]. It may be verified that our PBs reduce 
to those in [13] in the scalar case p = s = 1 upon constraining to w = 0. 

Remark 4-3. Notice that for r = p = s = 1 the AKS hierarchy on M. is the gl2 version 
of the well-known AKNS hierarchy. For r = 1 and arbitrary p, s (see Remark 2.2.), it 
is reasonable to call the system on M. a generalized AKNS hierarchy. For the generalized 
AKNS hierarchy the pseudo-differential Lax operator associated to C E M becomes just the 
operator AK (4.22). For this reason we can call AK in the factorization (4.21) the 'AKNS 
factor". (It is an easy exercise to directly verify the equivalence between the respective 
formulae (1.2-3) and (4.34-35) for r = 1.) In the simplest case r = p = s = 1 the connection 
between the AKNS hierarchy on M and the constrained KP hierarchy on M was observed 
in [12,13,25] too. The nonlinear Schrodinger (NLS) hierarchy results from constraining the 
AKNS hierarchy, and it has many generalizations [26] . The connection between generalized 
AKNS (iV-wave) and NLS systems and constrained (matrix) KP systems given by the above 
result and the results in Section 5 can be extended to more general cases than those treated 
in this paper. 
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To summarize, in this section we investigated the properties of the mapping 7r 
(4.15) and found that the push forward of the bihamiltonian structure on M c /Af in- 
duced by the Drinfeld-Sokolov reduction is the Gelfand-Dickey bihamiltonian structure 
on M = tt (M c /M). We presented the explicit formulae of the PBs on M c /Af. We have 
also established the factorization (4.21) of the Lax operator L G M using the Miura map 
\x : 6 — > .Mds mentioned in the paragraph following (2.24). This Miura map together 
with (2.24), and formula (4.34) of Theorem 4.4 provide alternative means for computing 
the reduced second PB. The commuting Hamiltonians furnished by the local monodromy 
invariants of C G M. c define generalized KdV and modified KdV systems in terms of the 
gauge slices .Mds and 6, respectively. 

5. Local monodromy invariants and residues of fractional powers 

In Section 4 we established a relationship between the Poisson brackets on the reduced 
phase space M re d — M. c /M and the Gelfand-Dickey Poisson brackets onMcA Our next 
task is to characterize the Hamiltonians generated by the local monodromy invariants of 
C G Ai c - These Hamiltonians, which define the commuting hierarchy of evolution equations 
on .M r ed, turn out to admit a description purely in terms of the Lax operator L G M 
attached to C by the elimination procedure, L = tt(C). Namely, the Hamiltonians defined by 
the local monodromy invariants of C will be identified in terms of integrals of componentwise 
residues of fractional powers of the diagonal PDO L obtained by diagonalizing L in the PDO 
algebra A. This identification will result from computing the local monodromy invariants 
of C G .Mds m two alternative ways: first using the procedure of (1.10-11) outlined in 
Section 1 and then using the diagonalization of L combined with a reverse of the elimination 
procedure. The same method was used in [1] , but the presence of the singlets in the partition 
(2.1) gives rise to complications requiring a non-trivial refinement of the argument. 

5.1. Local monodromy invariants and solutions of exponential type 

We wish to compute the local invariants of the monodromy matrix T associated to the 
linear problem 

C^ = (d x +j(x)+A)^(x) = 0, (5.1) 

where j(x + 2n) = j(x) since C = (d + j + A) G M c . If ^> : (R — > GL n is a solution of 
(5.1), which means that the columns of the matrix \& are a complete set of solutions, then 
the monodromy matrix is given by 

T = ^(2tt)*- 1 (0). (5.2) 
Following the procedure outlined in Section 1, perform the transformation 

C^C = e E Ce~ E = (d + h + A), * = e s ^, (5.3) 
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where 

H G (Im(ad A)) <0 , h G (Ker(adA)) <1 , (5.4) 

and the subscripts refer to the grading d in (2.7). The fact that S and h are uniquely deter- 
mined differential polynomial expressions in the components of j implies that S(j'(x)) and 
h(j(x)) are periodic functions of x G IR . It follows that the invariants of the monodromy 
matrix T that we are interested in are the same as the invariants of the transformed mon- 
odromy matrix 

f := ^(2tt)^- 1 (0) = G , - 1 (0)TG(0), (5.5) 

with the definition 

G(x) :=exp(-S(j(x))). (5.6) 



Using the notation (2.2-3), h in (5.3) may be written as 
where A(j) and D(j) are uniquely determined series of the form 

oo p oo 

A(j) = J2J2 h kAJ)K k D(j) = ^A- fc e r+ i, r+ i®D fc (j), (5.7b) 



fc=0 i=l k=0 



with hk,i(j(x)) G (C and D}~(j(x)) G gl s . A basis of the centre of Ker(adA) C i(gl n ) is 
given in (2.5) and (2.11). The Hamiltonians defined by the procedure in (1.10-11) are then 

Hk,i{j) '= dxh k!i (j(x)), i = l,...,p, k = 0,1,2,... (5.8) 
Jo 

and 

/■2tt 

E k (j) := / dxtrD k (j{x)). (5.9) 
Jo 

Taking the trace of eq. (5.3) using the identity A~ kr = X~ k l r , we obtain the equality 

r-2-K V 

E k (j) = 4,0 / dxtrj(x) -rVfffc riJ (j). (5.10) 
J ° i=l 



Since j i— > dxtrj(x) defines a Casimir function with respect to both Poisson brackets 
on M re d, this equality means that the complete set of independent Hamiltonians associated 
to the centre of Ker(ad A) is given by the H k ,i(j) above. General arguments that go back 
to the r-matrix (AKS) construction (see e.g. [10]) guarantee that the Hamiltonians H kyi 
are in involution (commute among themselves) since they can be interpreted as particular 
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monodromy invariants. To explain this interpretation, notice that the transformed linear 
problem = has the solution 

where 

-xA r 0T-^^2 d£h Ki {j{£))k- k ®e i A , (5.12) 

— i — 1 / 

and 

- y o dej]A-*i? fc o'(0)J- (5-i3) 

The corresponding monodromy matrix is 

T=f^ f A. (5.14) 
V * 22 (27r) J y 1 

To diagonalize T, consider an rth root ( of A, £ r = A, and define 

( a := (uj a with oo := exp {2m /r) . (5.15) 

The matrix A r is conjugate to 

A r :=diag (Ci,C2,-..,Cr). (5-16) 

In fact, we have 

A r = SA r S-\ (5.17a) 

with 

Sa b =^(C b r\ (S^a^^ita) 1 -*. (5.17b) 

Using this conjugation the upper block ^n(27r) of T becomes diagonal, 

(oo V \ 
-27rA r ® T - Hk >ti) ~ K r k ® e M i S ~ l ® ! P ) • (5-18) 

k=0 i=l / 

Hence, up to the constant — 2-7rA r <g> T, the Hamiltonians Hk,i(j) can be identified as ex- 
pansion coefficients defining the expansions of logarithms of certain eigenvalues of the mon- 
odromy matrix around ( ~ oo. As usual, this expansion has to be interpreted as an 
asymptotic — or formal — series. It is clear from (5.13) that the spectral invariants deter- 
mined by the "lower block" ^22(2^) of the monodromy matrix are in general, except the 
functionals Ek(j) given above, non-local functionals of j. 



28 



We have seen that the local monodromy invariants H k;i (j) associated to the centre of 
Ker(adA) are determined purely in terms of the upper block ^n(27r) of the transformed 
monodromy matrix T. To see what this means in terms of the original linear problem (5.1) 
consider the solution \& given by 



:=G(x)V(x) 



S 




1- 



(5.19) 



Using a block notation similar to (2.2-3), 

*w=(I"S !! 2 i:!V cm 



-s(i(*)) 



'Gu(x) G 12 (x) 
G 21 (x) G 22 (x) 



we have 



and 



*21 



*12 
^22 



G 2 i^iiS 



Gi 2 ^ 22 
G 22 fy 22 



(5.20) 



(5.21a) 



(5.21b) 



It follows from (5.4) that the infinite series S(j(x)) contains only non-positive powers of A 
and therefore the entries of the matrix G are given by similar series. This together with 
(5.12-13), (5.17) implies that the column vector solutions of the linear problem £ip = 
comprising the matrices in (5.21a) and in (5.21b) are qualitatively different. They are 
different in the sense that — apart from the constant S(Q <8> 1 P that we included in the 
definition of \& for later convenience — the columns of the matrix in (5.21a) have the form 
of descending series in non-positive powers of £ multiplied by a leading term of the type 
g-^CbT w hii e the columns of the matrix in (5.21b) do not contain such a leading term only 
a descending series in non-positive powers of (. We refer to the series solutions containing 
a leading term e~ x< * hT as "solutions of exponential type" . The solutions of exponential type 
contain all information about the local monodromy invariants since the matrix ^>h{2ti) in 
(5.18), whose eigenvalues generate the Hamiltonians H^ ^j), is conjugate to the matrix 



* n (27r) (^n(O))- 1 = Gn(0)*ii(27r) (Gii(O))" 1 



(5.22) 



It is convenient to write the (rp + s) x (rp) matrix in (5.21a) in the following detailed form: 



*21 



^ ••• V 2 



VV 1 Vr 
Xcj) 1 ■■■ (p r ) 



(V>a e gl p , (p b e mat(s x p)). 



(5.23) 
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Putting the above formulae together we have 



C) = ^ E A ) (^) c_1 ex P Cb) - *0T) , 

r 7 1 ( 5 - 24a ) 
() = 1=J2 G 2i fa A) (C6) c_1 exp (-V(x, C b ) - x( b T) , 

* c=l 

where and G^i are p x p and s x j) matrices, respectively, and £>(£, C) is the p x p 
diagonal matrix series 

D(x,C) = ^r fc / ^diag(Vi(j(0),---,^, P (j(0))- (5.24b) 
fc=o - 70 

Observe that all the gauge invariant components ip\(x,C) can be obtained from ip[(x,() 
simply replacing the argument £ = ( r by £5. By the reverse of the elimination procedure (see 
(4.11)), this means that the complete set of solutions of exponential type can be recovered 
from ipi(x,£). In particular, the local monodromy invariants H ki in (5.8) can be read off 
from the relation 

00 

i; r i(x + 2tt, C) = C) exp(-27rCr - ]T C'diag (H k , u . . . , H k , p )) , (5.25) 

fc=0 

which is a consequence of (5.24). This relation will play a crucial role in finding the link 
between the local monodromy invariants and the residues of fractional powers, which is the 
ultimate aim of the present section. 

In this subsection we obtained the matrix solution (5.23) of exponential type to the 
linear problem (5.1) by transforming £ to £ and imposing on the solution \I/ (5.11) of = 
the condition ^ = l n at x = 0, see (5.12-13). We noticed that the matrix solution in (5.23) 
is determined completely by the block We then observed that the block ipl(x, () directly 
encodes the Hamiltonians of our interest, the Hamiltonians H k)i given in (5.8), through the 
relation (5.25). In the next subsection we will consider the consequence of looking at ip{ 
from a different point of view, namely as a solution of equation (4.12), which was obtained 
via the elimination procedure. The result will be an explicit connection between the family 
{Hk t i} of local monodromy invariants and the components of the residues of the fractional 
powers of the diagonalized form of L. 
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5.2. Solutions of exponential type and residues of fractional powers 



Consider the p x p matrix PDO 

L = A r d r + u 1 d r ~ 1 + u 2 d r ~ 2 + --- + u r - 1 d + u r + z+(l s d + w)- 1 z-, A = -r _1 , (5.26) 

attached to £ E M c by the elimination procedure. The strategy of this subsection will be 
to determine ip{ in (5.24) as a solution of the linear problem 

Lipi = \ip! (5.27a) 

given by a (asymptotic or formal) series of the form 

M*X)= [f^Xk(x)C k )e- x(:r with XkeC°°(R,gl p ), det ( Xo (x)) ^ 0. (5.27b) 
\fc=o / 

The elimination procedure implies that ipl(x, Q in (5.24) is a solution of (5.27a). We shall 
see below that ip[ (x, Q can be expanded in the form given in (5.27b) and that the solution 
of (5.27a), (5.27b) is essentially unique. 

That iftl in (5.24) can be expanded in the form given in (5.27b) can be seen by in- 
spection. The key step is to check that the series S defining G = e~~ = l n — S + ^S 2 • • • 
in (5.20) contains only negative powers of A in its first row due to the grading condition 
(5.4), which implies that the leading term of comes from the unit matrix contained in 
the c = 1 contribution in the sum in (5.24a). Computing the first term of the "abelianised 
current" h(j) in (5.7), one obtains 

^,i(j) = ^((-r) r «i)„, (5.28a) 

where u\ is the gauge invariant component of j entering the Lax operator L (5.26) attached 
to C = (d+j + A). It follows that when rewritten as a series of the form (5.27b), the leading 
term xo of V'i m (5-24) is given by 

xo(x) = exp (-±(-ry J* <% («i(0)dtag) > ( 5 - 28b ) 

and indeed has non-zero determinant. Incidentally, the constant factor S(() was inserted in 
the definition (5.19) to set the leading power of £ in (5.27b), which multiplies the product 
of xo (5.28b) and e" xCr , to be C°- 

Below our aim is to determine ipi(xX) from equations (5.27a) and (5.27b). To make 
precise the meaning of equation (5.27a), which has been derived by formally applying the 
elimination procedure, we note that an arbitrary pxp matrix PDO a = ^ 4 ai(x)d l acts on 
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a series of the form ip±(x, () in (5.27b) as follows. Defining the action of d l for any integer 
% on e" xCr by (c^e" 3 ^) := (-Cr) J e- xCr one first associates the PDO x to tpi (5.27b) by 
writing 

oo 

Mx,() = (xe"< r ), i.e. X (x,d) = J2xk(x)A- k d- k . (5.29a) 

£;=0 

Then (aipi) (x, () := (/3e _x< > r ) with f3 = a\ being defined by the composition rule of PDOs 
(as in (4.2)). To avoid confusion, we stress that here the coefficienst of the PDOs %, a, f3 
are not required to be periodic functions on IR . When equation (5.27a) is understood in 
this sense it is easily seen to be equivalent to the "dressing equation" 

X~ X L X = A r <9 r . (5.29b) 

This reformulation of (5.27) is well-known, see [27] and references therein. Using the refor- 
mulation (5.29) and the fact that A r has distinct, non-zero eigenvalues, it is not hard to 
verify 2 that (5.27a) uniquely determines the series solution ipi(x, Q of the form (5.27b) up 
to multiplication on the right by a constant (^-independent) series c(Q of the form 

oo 

c(C)=J2 c ^~ k with det(c o )^0, (5.30) 

k=0 

where all the Ck are p x p diagonal matrices. 

Thanks to the above uniqueness property, the following procedure may be used to 
determine the series tpi(x, (). First we determine a p x p matrix PDO g of the form 

oo 

g = l p + J2^d-\ (5.31) 

i=i 

with periodic coefficients, gi(x + 2tt) = gi(x), such that 

L = gLg- 1 (5.32) 

where L is diagonal, i.e., 

oo 

L = A r d r + ^d r ~\ a, : all diagonal. (5.33) 

i=i 

Since A r is a diagonal matrix with distinct, non-zero entries, then if we require the g^s to be 
off-diagonal matrices we can recursively determine both the fi^'s and the a^s by comparing 

2 A proof can be found in [28] for the case when (wi)diag = and Xo = lp> but these 
assumptions can be dropped. 
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the two sides of Lg = gL term-by-term, according to powers of d. The solution is given by 
unique differential polynomial expressions in the coefficients defining the expansion of L in 
d. For instance, we have 

Al = Wdiag- ( 5 - 34 ) 

Then we consider a p x p matrix series 

V>i(z,C)= (jtc k Xk(x))e- xCr with det(xoOr))^0, (5-35) 
\k=o J 

which satisfies the equation 

(L^ 1 )(x,C) = A^ 1 (x,C), (A = C r )- (5-36) 
At the end of the procedure, we find the desired solution of (5.27) from 

Mx,()--= (gi>i)(x,0- (5-37) 

We know that up to a diagonal constant matrix c(Q of the form (5.30) the series solution of 
(5.27) determined by this procedure must coincide with ipl(x,£) given by equation (5.24). 
To proceed further we need the following identity. 

Proposition 5.1. The operator l p d can be expressed as 

00 \ - fc 

l p d = -TL 1 ^ + FkA k (L 1 ^ J , (5.38) 

fc=0 

with l}l r being defined by 

oo 

(L 1 /r)r = £ ) L 1 ^ = A8 + J2 b * d ~ l for some h, (5.39) 

i=0 

and uniquely determined p x p diagonal matrix valued functions Th, which satisfy 

F = -^(-r) r ai and J dx (kF k + (-T) k res(L fc / r )) (x) = for k > 0. (5.40) 

The above proposition, which was crucial in [1] for obtaining results analogous to those 
under consideration here, is taken from [4] and is originally due to Cherednik [29] (see also 
[27,30]). Since tfji is uniquely determined by (5.35-36) up to multiplication by a diagonal 
constant matrix c(() of the form given in (5.30), equation (5.36) implies 

(L^iXz, C) = Cl?i(*,0- (5-41) 
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This together with (5.38) leads to 

oo 



fc=0 

oo 



^Vi(^C)=exp(-<r + ^A fc r fc / ^^(0)^1(0,0, (5-42) 

fc=0 >/o 

where V>i(0, £) is an arbitrary diagonal matrix series of the form given in (5.30). Combining 
this with (5.34) and (5.40) results in the crucial relation 

Mx + 2n,C)=M^0r(C) (5.43a) 

with 

-2<r - ±-L y o dx ( Ul (x)) di&g - E V i, dxres r /r ) {x) ) ■ 

(5.43b) 

Observe now that ip\(x,Q defined in (5.37) satisfies 

^i(x + 27r,C) = ^i(^CMC) (5-44) 

with the same t(£), because the coefficients gi(x) defining the PDO g are periodic functions 
of x. This immediately leads to the following result. 

Theorem 5.2. The set of commuting Hamiltonians provided by the local monodromy in- 
variants Hk,i(j) given in (5.8) is exhausted by the Hamiltonians Ho^j) together with the 
Hamiltonians defined by the residues of componentwise fractional powers of the PDO L 
(5.33) obtained by diagonalizing the Lax operator L (5.26) attached to C = (d + j + A) by 
the elimination procedure. More precisely, 

diag (#0,1, • • • , H , p ) = / dx («i(x)) diag , 

J° ' ( 5 - 45 ) 

f 1 ( ~ \ 

diag(H kjl ,...,H kjP ) = -J dxres[L k/r ) (x) , for k > 0. 

Proof. The statement follows by comparing (5.25) with (5.44) taking into account that 
ipi(x, C) = ipi(x, C) C (C) where c(£) is a diagonal constant matrix of the form (5.30). Q.E.D. 

Finally, we shall write down the evolution equations generated by the above Hamilto- 
nians on A4 Te d = Ai c /Af. For this it is convenient to consider an arbitrary p x p diagonal 
constant matrix Q and associate to it the Hamiltonian 

v 

H k ■= rY,Qa H Ki- (5-46) 
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Observe that has the form (4.39) since 



HQ(l,z + ,z-,w)=H%(L) 



Qi 



with 



tQi 



2n 



dxtr((-T) r Q Ul (x)) , 
k > 1. 

Using (4.40) the gradients of these functions are found to be 



H^(L) = -Tr[QL k / 1 



(5.47a) 



(5.47b) 



XT/® 

■f- = (-ry Q«r r , 



5L v ' ^ 5L 
where g is the PDO (5.31) which diagonalizes L and 



<)V i- = gQg-'L^/r, k > 1, 



L l/r 



iL l ^g-\ 



VZ. 



(5.48) 



(5.49) 



Note that L l l r commutes with gQg~ x . We let X£, q denote the Hamiltonian vector field 
associated to by means of the first and second PBs on A4 re d given in Theorem 4.4. 
Inserting the gradients (5.48) into formulae (4.43) and (4.44) we obtain the following result. 

Corollary 5.3. The Hamiltonian vector fields Xji. q defining the local hierarchy of compat- 
ible evolution equations on M. c /N take the form 

n, Q (L) = n +r , Q (L) = [(gQg- 1 L k / r ) + ,L , 

n, Q (z+) = Xl + r, Q (z+) = Po (gQg-'L^z+W- 1 ) W, (5 _ 5Q) 

X 1,q( w ) = X k+r,Q( w ) = °> Vfc = 0,l,.... 

In particular, the flows are bihamiltonian. In terms of the notation 



K 

{gQg- 1 L k ' r ) + :=Y,^ 



(5.51) 



i=0 



and the covariant derivatives defined in (4- 38b) the second and third equations in (5.50) can 
be rewritten as 



i=0 

k 



(5.52) 



i=0 
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This completes our general analysis of the KdV type hierarchy resulting from the gener- 
alized Drinfeld-Sokolov reduction defined in Section 2. The evolution equations associated 
to the vector fields in (5.50) define natural "covariantized matrix generalizations" of the 
constrained KP hierarchy considered previously in the literature (see [12,13,14] and refer- 
ences therein) in the scalar case p = s = 1 with the constraint w = 0. Notice that w and the 
diagonal components of u\ (the subleading term of L in (5.26)) do not evolve with respect 
to the flows determined by the vector fields in (5.50). The reason for this is the fact that 
(ui)a for % = 1, . . . ,p and the components of w are Casimir functions with respect to the 
reduced first PB, and generators of residual symmetries with respect to the reduced second 
PB. Indeed, evaluating the gauge invariant current component $i given by (3.3a) in the DS 
gauge (2.19), we see that <E>i is proportional with {u\)a. 

Remark 5.1. In the above we focused on the local mono dromy invariants, but it is clear that 
there exist also non-local monodromy invariants, generated by the eigenvalues of the lower 
block ^22(2^) of T in (5.13-14). The possible role of these non-local monodromy invariants 
concerning the integrability of the system should be further studied. In particular, it is not 
obvious to us whether or not these are to be regarded as independent from the commuting 
local Hamiltonians that we described. (An argument valid in finite dimensions would imply 
that the gradients of all invariant functions on the dual of a Lie algebra belong to the 
centre of the centralizer.) It would be also interesting to know whether or not, apart from 
functions of the current components $i and w in (3.3a) and (3.3b), there exist further (local) 
commuting Hamiltonians that are not generated by monodromy invariants. 

6. Examples for r = 2 and factorizations of the AKNS factor 

The purpose of this section is to elaborate the simplest examples for r = 2 in order 
to illustrate the preceding results and to make the abstract developments more concrete. 
A particularly interesting result that will be obtained from considering the r = 2 case is a 
further factorization of the AKNS factor AK (4.22) appearing in the modified KdV Lax 
operator (4.21). In fact, we shall derive the factorization 

K = (l p d + a - 6(1,0 + cO^c) = (l p d + 0_i)(l p - A/3(l s d + tf + l^P)' 1 !), (6.1) 

and a similar one containing a factor (l p d + $1) on the right, in terms of the new variables 

j3 G mat(p x s), 7 G mat(s x p), di G gl , i = ±1, $0 G gl s . (6-2) 

Substituting this factorization of K in (4.21) leads to a second modification of the gener- 
alized KdV hierarchy for any r. The variables in (6.2) will be seen to parametrize a new 
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gauge, the so called G'-gauge, which will arise as a gauge section in an alternative version of 
the Drinfeld-Sokolov reduction leading to the same reduced phase space M c /Af (2.16-17). 
In the first part of this section, we shall obtain the explicit form of the reduced second PB 
in terms of the G'-gauge as well as the G-gauge, the G-gauge and the DS gauge defined 
in Section 2, and shall also present the local Poisson maps given by gauge transformations 
between these gauges. The simplest evolution equations of the generalized (modified) KdV 
hierarchies for r — 2 will be worked out in the second half of the section. 

6.1. The reduced second PB, Miura maps and factorizations 

Below we shall use 3 different gradings of Q := gl n , n = 2p + s, defined by the adjoint 
actions of the following 3 matrices: 

I = ^diag(l p ,-lp,0 8 ), H = diag (l p , P , S ) , H = diag (0 P , -l p , S ) . (6.3) 
For example, we have 

G = Q\ + 0i° 1/2 + Go" + G{° /2 + G{° , (6.4) 

where G^° is the eigensubspace of ad/o with eigenvalue k. Using a 3 x 3 block-matrix 
notation, the grades assigned to the entries of m G Q are given by 

/ 1 1/2 \ / 1 1\ / 1 \ 

m h) =1-1 -1/2 , m H = -1 , m H = -1 -1 , (6.5) 
\-l/2 1/2 0/ \-l 0/ \ 1 / 

in correspondence with 7 , H, H above. Previously we made use of the Ji"-grading to 
define the generalized Drinfeld-Sokolov reduction. Now we define an alternative version of 
it using the Icrg ra ding. In this version of the reduction procedure we first introduce the 
new constrained manifold M.' c C M. given by 

M' c := { C = d + j + A I j : S 1 - G% /2 }, (6.6a) 

where the explicit form of A and j is 



(6.6b) 



Then we factor ize M' c by the new "gauge group" Af' generated by G< - An element g e Af' 
acts on M' c via the gauge transformation 

/l p \ 

g:C^gCg-\ g=\* l p * : S 1 ^ exp (^° ) . (6.7) 
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Comparing the -ff-grading with the / -g r ading one sees that M c C M' c and M C A/"', where 
A/f c and A/" are given in (2.16-17). This enlargement of the constrained manifold and the 
gauge group has been defined in such a way that the factor space is unchanged, 

M re d := M c /M ~ M'JW. (6.8) 

This follows since there exist unique gauge transformations that bring an arbitrary element 
of the constrained manifold to the DS gauge A/Ids, 

/ \ 

Ads = d + \v 2 v 1 C+ + A for £ DS E M DS , (6.9) 
VC- wj 

in either of the two reduction procedures, as one may verify by computing the (differential 
polynomial) formulae of these gauge transformations. 

The interest of the new reduction procedure is that it admits the gauge section ©' C M' c 
whose general element Cq> E 0' takes the form 

M o p\ 

£ e , = d+ 0_i + A, (6.10) 
V 7 O / 

that is, ©' is defined by restricting j to take its value in {Qq° + Q[° 2 )- We a ^ so nave the 
©-gauge and the ©-gauge defined by the sub manifolds C Ai' c , © C A4' c where the 
constrained current j is restricted to Qq and Qq , respectively. We may parametrize the 
general elements Cq E and £@ E as follows: 

#i 0\ fa b\ 

C e = d+\ a 6 +A and £q = 5 + °-i + A. (6.11) 
c dj \c d J 

Let now F, H be gauge invariant functions on and consider the function 

P:={F,fT} 2 , (6.12) 

determined with the aid of the usual extension-computation-restriction algorithm, which is 
also a gauge invariant function on M.' c . Of course, the functions F, H, P can be recovered 
from their restrictions to any of the above mentioned gauge slices, denoted respectively as 
Fq>, Fq, Fq, Fds and similarly for H and P. The function Pq> may be determined in 
terms of Fq/, Hq> using the following simple formula: 



P ,(0_i, 0o,0i, P,l)= V / tr (0, 

^=-i,o,i Jsl V 




6F& SH e 
5F e , 5H e 

~Jy W <5 7 5/3 



(6.13) 
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The formula on the right hand side describes the reduced second PB in the variables 
parametrizing ©'. Here the functional derivatives are defined in the standard way using 
the scalar product provided by the ordinary matrix trace. The derivation of this formula is 
rather straightforward and we omit it. The formula of the reduced second PB in terms of 
the variables parametrizing the 6-gauge can be read off from (2.24), 



i=0,l ' 



SF e SHq 
86, ' 86 i 



SFq (5H e y 
86, V S9i J 



6 := 



and in terms of the O-gauge the formula is similar, 



i=0,-l 



tr Oi 



SFq 8H @ 
66, ' 86, 



SFq (8H t 



86, V 86, 



6o 



a b 
c d 



a b 
c d 



(6.14) 



(6.15) 



Now we wish to present the PB in the DS gauge. It is convenient to do so in terms of the 
variables u,v,y±,w that appear in the Lax operator 



L = A 2 (d 2 + ud + v + y+(d + w^yJ) , 
attached to £ds (6.9) by the elimination procedure, i.e., 

u = -ruiA, v = -Tv 2 , y+ = T( + , y_=C-, 



(6.16a) 



(6.16b) 



since then the entries of the constant matrix A = — r _1 will not appear explicitly in the PB. 
To simplify the formula, which can be obtained by a direct computation or by specifying 
(4.34), we introduce the notation 



8H DS 



8H DS 



8H DS 



8H DS 



8H DS 



(6.17a) 



bu bv by+ by- bw 

The function Pds(u, v, y±, w) defined by restricting P in (6.12) to the DS gauge is given by 



DS 



tr 



S 1 



bu 



(6.17b) 



where the formula of the Hamiltonian vector field X 2 ^ is found to be 



X.% (u) = [£„, u] - 2& + v] + C - (W + S-v- - 

X 2 ^) = v] - C - < + <vv - vi v u + < - - u^ v u)' + < + (M' 

+ C + y+v(Z+) - y+t+u + (C-y-Y + + <-y- + y+y-}, 

X. 2 H (y+) = £ u y + + u£ v y+ + V(^y+) + V 2 {t-) + uV{^_) + <_ - y + £ w , 
X.%{y_) = -y_£ u + V(y_£ v ) + y_g v - y_^ v u - V 2 (£ + ) + V(£+u) - £+v + 
X 2 H (w) = i + y + - y-i- + [£ w , w] - C, (6.17c) 



39 



with the covariant derivatives defined in (4.38b). 

Using (6.7) and (6.9-11), a straightforward inspection shows that there exist local gauge 
transformations (for which the entries of g G M' acting on£ = <9 + J+ A are differential 
polynomials in the components of j) implementing the following changes of gauge: 

i) The mapping v : Q' — > ©, where Cq = v(jCq') = g v C,®ig~ x is given explicitly as 

0i=0i, a = tf_i-A/3 7 , d = tf + 7A/3, b = A/M - Af3<yA[3 + A/3', c = 7 . 

(6.18) 

ii) Analogously, the mapping v : 0' — > 0, where Cq = u(Cq') = g^C^g^ 1 is given by 

0_i = <0_i, a = #i-/?7A, J=^o + 7A/3, 6 = /?, c = 7Ai?i -tfo7A-7A/37A-7'A. 

(6.19) 

iii) The mapping /i : — > -Mds 5 where £ds = A*(^g) = Q^eO^ 1 is defined by 

vi = a-A0ir, w 2 =aA0i -6cA + A0i, C+ = &, C- = cA0i - cfcA - c A, u; = d. 

(6.20) 

iv) Similarly, the mapping j2 : © — > .MdSi where £ds = A*(£e) = 9p, J ^e9ji 1 rea ds 

vi = 0_i-Aar, u 2 = 0_i Aa- Afec+Aa', C+ = 0-iA6- AM + Aft', C- = c, w = d. 

(6.21) 

The construction ensures that these "generalized Miura maps" are in fact Poisson maps 
with respect to the PBs given above that realize the reduced second PB in terms of the 
respective gauges. In accordance with what is expected from Miura maps, the inverses of 
these maps are not single valued and are non-local. Incidentally, the map v : ©' — > 
(similarly the map v) provides a realization of the gl p + s current algebra appearing as a 
term in the PB (6.14) on in terms of the gl p © gl s current algebra and the "symplectic 
bosons" /?, 7 appearing in the PB (6.13) on ©'. This realization of the gl p + s current algebra 
is analogous to the well-known Wakimoto realization of the gh current algebra (see [31] for 
an elaboration of this observation) . 

Now we come to the relationship between the above Miura maps and the factorizations 
of the Lax operator. Observe first that the linear problem Cip = has similar properties 
for C G Ai' c as for C G M. c - It is covariant under 

g : C I— > Ce~^ ', ip i— > gip for any g G A/ 7 , (6.22) 

and the component ipi of ip = is invariant. Thus we may apply the elimination 

procedure and on account of (6.8) the resulting image of M.'JM' will turn out to be the 
same manifold of PDOs which was obtained in Sect. 4. At the end of the chain of Miura 
maps ©'-^©-^A^ds is the DS gauge to which the elimination procedure associates the PDO 
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L (6.16). We have also seen in Sect. 4 that applying the elimination in the 6-gauge leads 
to a factorization of this Lax operator, 

L = Lq = AKA(d + 9i) with K = (l p d + a — b(l s d + rf) _1 c) . (6.23) 

We now wish to show that performing the elimination procedure in the O'-gauge yields a 
further factorization according to 

L = Lq = L@> = A(l p d + 0_i)(l p - A/3(l s d + $ + 7 A/3)- 1 7 )A(l p a + (6.24) 

which amounts to the factorization (6.1) of the AKNS factor AK since 9± = $i by eq. (6.18). 
For Cq> E ©' (6.10) the linear problem takes the following form: 

(ipd + h) r p \ 

Ar (l p + 0_i) V2 =0. (6.25) 

7 (l a d + do)J \(f> J 

More explicitly, 

(i p a + ^ 1 )^ 1 + ^ + r^ 2 = o, 
(i p a + ^_i)^ 2 + Ar^i = o, (6.26) 

(1,0 + tf o )0 + 7^2 = 0. 
In order to obtain the desired equation, 

Le'Vi = A^i, (6-27) 

from the middle equation in (6.26) we must express ^2 in terms of •i/'i using the other two 
equations. To do this, we first "formally integrate" the last equation in (6.26) to yield 

<j> = -(l s c> + tfo)"V/>2. (6-28) 
Substituting this back, the first equation in (6.26) becomes 

(l p d + tfityi + (r - (3(l s d + ^o) _1 7)^2 = 0. (6.29) 

Finally, formally solving (6.29) for ip 2 , and then using the middle relation in (6.26), we 
obtain (6.27) with 

L e , = (A(l p d + #-!)) (l p + A/?(l s a + ^o) _1 7) _1 (A(l p d + #!)). (6.30) 
This is the same as L®> in (6.24) on account of the PDO identity 

(l p + A(3(d + ^o)" 1 ^" 1 = (1 P - Ap(l s d + tf + 7A/3)" 1 7 ) . (6.31) 
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The factorization (6.24) now follows since the covariance properties of the linear problem 
imply L = L® = Lq> for the Lax operators associated to gauge equivalent points of A4' c . 

Considering the alternative chain of Miura maps G'-^G-^.Mds m a similar manner, 
we obtain 

L = L e = L e , where L@ = A(l p d + 0_i)A (l p d + a — b(l s d + (J) _1 c) . (6.32) 

This is equivalent to the alternative factorization of the AKNS factor 

A (l p d + a - b(l s d + d)- 1 ^ = (l p - A(3(l s d + + 7A/3) _1 7) A (l p d + #i) . (6.33) 

It is worth noting that all the factors in equations (6.1) and (6.33) correspond to Poisson 
submanifolds in the PDO space A endowed with the quadratic Gelfand-Dickey PB (4.5). 



We here work out the first few Hamiltonians and evolution equations for the generalized 
(modified) KdV hierarchies with r = 2. For simplicity, we now restrict ourselves to the scalar 
case p = 1, where A can be set equal to 1. To compute the commuting Hamiltonians we 
need the expansion of L in powers of d, and from (6.16) using (4.33) we have 



6.2. Some explicit examples 



oo 



L = d 2 +ud + v + J2 L kd~ k with L k = (-l) k - 1 y + V k - 1 (y_). 



(6.34) 



fc=i 



The commuting Hamiltonians may be obtained from the residues of powers of L 1 / 2 i 
accordance with (5.47), where in the scalar case L = L and Q = 1. If we write 



m 




(6.35) 



then for the first three terms we find 




h = -L 2 --l 2 u+-u'l 1 --£ 2 1 --£> 2 . (6.36) 



Substituting the above formulae we get 




(6.37b) 



(6.37a) 



(6.37c) 
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Hs(L) = lj^es (L 3 / 2 ) = | J QuLi + L 2 + < + ^ + W £ 2 + ^ (6.37d) 

= y Qwy+y- - y+ v (y-) + \™>' + ^K) 2 + \v 2 - \u 2 v + ^u 4 ^j , 

Ha{L) = 1 y res (L 2 ) = y (L 3 + «L 2 + wL^ + vL x ) (6.37e) 

= y (y+^ 2 (y-) - uy+V(y_) + u(y+y-)' + vy+yJ) . 

These formulae describe the first few Hamiltonians in terms of the generalized KdV fields 
u,v,y±,w associated to the DS gauge. The expressions for the Hamiltonians of the modified 
KdV systems corresponding to the gauges 6 and ©' can be found as the pull backs of the 
above Hamiltonians by means of the Miura maps n in (6.20) and v in (6.18). For example, 
in terms of the variables 6\,a,b,c,d, the first three Hamiltonians of the modified KdV 
hierarchy associated to the 6-gauge read 

/j*H = J (a + 9^, n*Hi = - y (be + ^(a - 6>i) 2 ^) , /j*H 2 = J b{c' + cd - c0i) . 

(6.38) 

In terms of the variables (3, 7, 0j (i = 0, ±1), the corresponding Hamiltonians of the second 
modification associated to the ©'-gauge are given by 

1/VW0 = y + - Pi) ' ( 6 - 39a ) 

= -\J ((#1 " 0-i) 2 + 2 (^i + 0-i " 2tfo) 07 + 4/3' 7 - 3/3 2 7 2 ) , (6.39b) 

u*n*H 2 = J (d-i/3 - $o/3 - if + 13') (V + 7 (0o - #1 + /*y)) • (6.39c) 

The expressions of the higher Hamiltonians soon become very long. 

As an illustration, let us give the evolution equations in the different gauges determined 
by the Hamiltonian functions Hi, n*TL\ and v* [i*Hi with respect to the appropriate PBs 
given in Sect. 6.1. In the DS gauge we find the flow 

w = u = o, v = v' - \i u " + uu ')i y+ = v(y+) + \ u y+> y- = v (v-) ~ \ u y-- ( 6 - 40 ) 

In the G-gauge and in the G'-gauge we have 
d = 0, a = -$! = \{a-9i)', b = ]-ab+]-b9 1 -db+b' , c = ~ac~6ic+dc+cf , (6.41) 

— ■ Z Z Z Z 

and respectively 

O = -{M, 0±i = ±^(0i - 0-i) + \(Pl)', (6-42) 
13 = + 0-i - 20 o - 3/3 7 ) +/3', 7 = -^7(0i + 0-i - 20 o - Wl) + 1' ■ 
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It can be checked that the Miura maps 6'-^6-^.Mds map these flows one to another 
accordingly. In the same way one can in principle compute the evolution equations for the 
other Hamiltonians. Since the formulae are rather complicated we only give the next two 
evolution equations in the DS gauge. The simplest way to obtain these is to combine the 
formulae in (5.50) with the identities (L) + = d 2 + ud + v and 

(£ 3/2 ) + = d 3 + \ud 2 + \{v + \u 2 + \u')d + \{v> + 2y + y. + uv) - ±u 3 + \u" . (6.43) 

The flow determined by the Hamiltonian H2 (6.37c) through the PB (6.17c) satisfies 

w = u = 0, v = 2(y + y-)\ y + = V 2 (y + ) +uV(y+) +vy +7 y- = -V 2 {yS) +V{uyJ) -vy— 

(6.44) 

For the flow determined by 7i 3 (6.37d) by means of the PB (6.17c) we get w = u = 
together with the equation 

,00 o 1 

/// I /II 11 \ /ll 11 I 2 l\ 

v = —v + —vv + -{y + y- - y+y_) — -{uv +u v + uu v + -u v ) 

Zi Zi 41 

+ ±((u 3 )" - 2u"" + 3u 3 u' - 2uu"') - 3(wy+y-)' + \{uy + y-)\ 

V+ = V 3 (y+) + \uV 2 {y + ) + ^vV(y+) + ^u"D(y+) + ^u 2 V(y+) 
3 3 3 1 1 

+ ^(y+y-)y+ + + -^uvy+ - j^ u3 y+ + g u "y+> 

y_ = V 3 {y_) - \v 2 {uyJ) + ^D{vy_) + -V{u'y_) + ^V{u 2 y.) 

- 7^(y+y-)y- - \ v 'y- - \ uv y- + ^ u ^y- ~ \ u " y — 

This generalization of the KdV equation reduces to one of the equations in [13] when w = 
and u = 0. 

In the matrix case p > 1 the Lax operator L has to be diagonalized according to (5.32) 
to compute the Hamiltonians from the residues of fractional powers of L. The formulae 
that can be obtained straightforwardly do not appear particularly enlightening to us and 
they are too long to justify presenting them. 

7. Discussion: discrete reductions and generalized KdV hierarchies 

In this paper we have derived a Gelfand-Dickey type PDO model of the hierarchy 
resulting from generalized Drinfeld-Sokolov reduction in the case of £{gl n ) with the partition 

p times s times 

(fT~^~r, !,...,!), n=pr + s, (7.1) 



(6.45) 



44 



using the grade 1 semisimple element A in (0.6). The reduced phase space turned out to 
be the space of quadruples (£, z+, Z-, w) where i is a p x p matrix r-KdV type operator 
which is coupled to the fields z±, w. The compatible PBs and the commuting Hamiltonians 
obtained from the Hamiltonian symmetry reduction are given by Theorem 4.4 and Theorem 
5.2., respectively, with the corresponding evolution equations being described in Corollary 
5.3. These results extend the results of [1] on the matrix KdV system for which s = 0. 

Of course the Drinfeld-Sokolov method works in a more general Lie algebraic context 
than the case of £(gl n ) and a semisimple A of minimal positive grade that we have considered. 
However, it is worth noting that some of the systems that would result from the Drinfeld- 
Sokolov construction applied to £{Q) for Q a classical Lie algebra can be also obtained from 
"discrete reduction" applied to a generalized KdV system associated to gl n . In the rest of 
this section we wish to briefly explain how this comes about. 

Let us consider the discrete symmetries of the system described by Theorems 4.4 and 
5.2. More precisely, let us look for symmetries given by some involutive map a on the 
reduced phase space -M re d = M c /Af = {(£, z + , z_, w)}, 



a : M re d -> M red: 



= id, 



(7.2) 



which leaves the PBs invariant, 



{/ O^/iO a}* = {/,/*}* o a, i = 1,2, 



(7.3) 



for arbitrary functions /, h on Ai re d- We will take the following ansatz for a. Let m G GL p 
and let q G GL S , i.e., m and q are constant, invertible, respectively pxp and s x s matrices. 
Define the map a m , q '■ -Mred -Mred by 



>m,q 



/ £ \ / mtm~ x \ 
z + 



—mz^q -1 
qz\m~ x 
—qw t q~ 1 



where ft is given by the standard adjoint operation on the PDO space A, 

r r 

ft = (-i) r A r d r + ^2{-iy- l d r - l u\ for i = A r d r + u *d r 



i=l 



1=1 



(7.4) 



(7.5) 



It is not hard to verify that a m , q is a Poisson map with respect to the PBs given by Theorem 
4.4 whenever it maps the phase space -M re d to itself, which is ensured by the condition 



mAm" 1 = (-l) r A r . 



(7.6a) 



45 



The involutivity of cr m ,g leads to the conditions 

m* = e m m, e m = ±1, q l = e q q, e q = ±1, with e m e q = -1. (7.6b) 

Notice that if e m = — 1 then p must be even and when e q = — 1 then s must be even. For 
any natural numbers a and 6 define the a x a and 26 x 26 matrices rj a and O26 by 

a b 2b 

Va = ^2 e i,a+l-ii ^26 = ^ ej ; 2b+l-i ~ ^ ^,26+1-1, (7.7a) 

and let f a denote an arbitrary a x a diagonal, invertible matrix subject to 

Va^aVa = -fa, (7.7b) 

which means that f a is anti-symmetric under transpose with respect to the anti-diagonal. 
We have the following solutions for a m ^ q . 

Type i): r = 2p even, Vp, m is diagonal and g is arbitrary with e q = — 1 (s = 21 even). 

Type ii): r = (2p + 1) odd, p = 2k even, A is such that r\ p Ar\ p = —A, m = f p O p and q is 
arbitrary with e = — 1 {s = 21 even). 

Type iii): r = (2p + 1) odd, p = 2k even, A is such that r) p Ar) p = —A, m = ^ p r] p and g is 
arbitrary with e q = +1 (Vs). 

Note that rj p Arj p = —A requires A to have the form A = diag (Ai, . . . , A&, — A/-, . . . , — Ai), 
where p = 2k and AJ" 7^ =t^J 7^ for i ^ j since A r must have distinct, non-zero eigenvalues 
(cf. eq. (2.9) for T = —A -1 ). The symmetries of type i), ii), iii) are of course also available 
in the case of the corresponding r-KdV systems, where s = and the phase space is simply 
the space of operators £ on which o m ^ q acts by £ t— > mftmr 1 . 

Given an involutive symmetry a = cr m g , one finds that 

a : L^mL^m- 1 for L = £ + z + (l s d + w)~ 1 Z-, (7.8) 

which implies that the commuting Hamiltonians given by Theorem 5.2 admit a basis con- 
sisting of invariant and anti-invariant (that change sign) linear combinations with respect 
to a. On account of (7.3), if H o a = H then the Hamiltonian vector fields are tangent 
to the fixed point set M.° ed C M re d of a. 

The flows of a "discrete reduced hierarchy" may therefore be defined by restricting 
the flows generated on A4 re d by the cr-invariant linear combinations of the Hamiltonians 
in Theorem 5.2 to the fixed point set M.° ed . These flows are bihamiltonian with respect 
to the restricted Hamiltonians and a naturally induced bihamiltonian structure on M.° ed . 
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The induced PBs on M° ed are defined by restricting the original PBs of functions of o- 
invariant linear combinations of the components of £, z+, Z-, w — which may be regarded 
as coordinates on M° ed — to M° ed . For fixed p, r, s and a given symmetry type i), ii) or 
iii) the various possible choices of m and q defining cr m ,<2 are equivalent from the point of 
view of the discrete reduction. In fact, the fixed point sets corresponding to two different 
choices are always related by a Poisson map of M re d given by 

mz+q~ 



I £ ^ I m£m~ L ^ 



Z- 



qz-tn 1 



(7.9) 



V qwq 1 / 

with some constant matrices rh G GL P and q G GL S . 

As explained in particular cases in [19], the above discrete reductions are actually 
induced by the reductions of gl n to a simple complex Lie algebra of B, C or D type. 
Correspondingly, many of the generalized KdV hierarchies that may be associated to certain 
conjugacy classes in the Weyl group W(£?) for Q a classical simple Lie algebra by generalized 
Drinfeld-Sokolov reduction (see [6,7,8,9]) can be also obtained by discrete reducing either 
the matrix r-KdV hierarchy or its extended version associated to gl n . These KdV type 
hierarchies are associated to a regular or a non-degenerate, in the sense that it admits 
DS gauge fixing, graded semisimple element A of minimal positive grade from a graded 
Heisenberg subalgebra of £(Q) by means of the construction described in Section 1. (The 
construction described in Section 1 is the special case of the construction described in 
[6,7] for which the "coarser grading" of £{Q) is chosen to be the homogeneous grading.) 
Since the graded Heisenberg subalgebras are classified [17] by the conjugacy classes [32] in 
W(£/), we can label these generalized KdV hierarchies by the respective conjugacy classes in 
W(Cy). The conjugacy classes that occur here can be parametrized (as in [32,19]) by certain 
"signed partitions". Using this notation, we find that the discrete reduction operates on 
the generalized KdV systems associated to gl n according to the following reduction rules: 

p times 21 times p times I times 

°A,n 2l ■ (27~^~2p, T~^~l) G W(gl 2 (p P+ i)) =>- (pT^Tp, e W(C PP+Z ), 

2k times 21 times k times ^ I times ^ 

0-Afi 2fe ,fi 2! : (frrTT7?,l,..T7l) G W(gl 2 ( kr +l)) =>• (rv*^~r, 1, . . . , 1) G W(C fcr+ ;), 

2k times 2H-ltimes fc times ^ 1 times 

<7Ar, 2k ,r,2i+i : 1, . . . , 1 ) G W (gl 2{kr+l) + 1 ) =^ {f^~~T, 1, . . . , 1) G W(B kr+l ), 

2k times 21 times k times I times 

(J A^ 2fe ^ 2( : (rv^T^, !,•••,!) e W(^ 2 (fcr+z)) =>- (f^~?, 1, . . . , 1) G W(D fcr+; ), 

(7.10) 

where / > is arbitrary, r = 2p + 1 is odd. By detailed inspection, this result has been 
established in [19] for / = in all cases. Since the case of arbitrary / can be treated in 
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a similar way, we omit the proof 3 . It is worth noting that the above reduction rules are 
also valid in the generalized AKNS case (see Remarks 2.2. and 4.3.) as well as for the 
non-abelian Toda systems corresponding to the generalized KdV systems. 

Finally, we wish to remark that not all generalized (non-principal) KdV hierarchies 
based on a classical Lie algebra are discrete reductions of hierarchies associated to gl n . 
This is apparent from the classification of graded regular elements [18,19], since certain 
"primitive" cases exist for Q = D 2p in correspondence with the conjugacy class (p,p) in 
W(D 2p ). The two negative cycles (p : p) also appear as building blocks of more complicated 
regular conjugacy classes in W(D n ) and in W(£? n ). A PDO description of the KdV systems 
associated to these conjugacy classes in the Weyl group by generalized Drinfeld-Sokolov 
reduction is not known. 
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3 The 1 = 1 case of the -Dfcr+Z series is exceptional since in this case the Drinfeld-Sokolov 
reduction applied to D^r+l gives more commuting Hamiltonians than the discrete reduction of 
the corresponding hierarchy based on g^kr+2- The relevant grade one semisimple element A G 
£(gl2kr+2) n£(-Djtr+l) is a regular element of £(Dk r +l) but is not a regular element of £(gl2kr+2) ■ 
The abelian subalgebra Ker(adA) C £(Dk r +i) contains extra generators in addition to those 
generators of the centre of the centralizer of A in £{gl2kr+2) that belong to £(Dk r +i) C £{gl2kr+2) ■ 
These yield the extra Hamiltonians by means of the Drinfeld-Sokolov construction. The form of 
these Hamiltonians in terms of the reduced phase space variables is not known to us. 
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Appendix A: The Poisson submanifold M K C A 



In this appendix we consider the mapping rj : Go — > C .4 given in (4.25) and show 
that Mx C A is a Poisson submanifold with respect to the quadratic Gelfand-Dickey PB 
and r\ is a Poisson mapping as stated by Proposition 4.1 in Sect. 4. 

It will be convenient to reformulate the statement of Proposition 4.1 in terms of Hamil- 
tonian vector fields. Let H be an arbitrary function on A, the corresponding Hamilto- 
nian vector field on A defined by means of the quadratic Gelfand-Dickey PB, 

x ^ = { l jz) + l - l {jz l ) + ' V£e * <A1 » 

and the restriction of X-^ to Mk C A. Then consider the function H on Go given by 

H = Ho V , (A.2) 

and denote by Y# the corresponding Hamiltonian vector field on Go defined by means of 
the current algebra PB (4.26). If 

J) (A3) 

denotes the arbitrary element of the space Go = gl p + s as in (2.22b), (2.23) and the gradient 

m / SH SH \ 



is defined in the usual way, one has 

Y H (8 ) = 



5H 



Clearly, the claim of Proposition 4.1 requires the vector field X-^ to be equal to the push 
forward of the vector field Y# by r), i.e., we must verify the equality 

Xh(v(0o)) = V*(Yh(0o)) V^oGGq. (A.6) 

In particular, the required equality implies that X^ is tangent to Mk (which is just the 
statement that Mk C A is a Poisson submanifold) and that i]*(Yh) gives a well-defined 
vector field on Mk- 

The formula for 77, 

K = v (6 ) = (d + a)-b(d + d)- 1 c, (A.7) 
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and (A. 5) permits us to write down the explicit form of the r.h.s. of (A. 6) as follows: 



V* (Yh (0 O )) = Y H (a) - Y H (b)(d + d)-^ - b(d + d^Ynic), 



(A.8) 



where we have 



Y H (a) 



SH 

5a 



,a 



(SH\' 5H b 5H 
\ 5a J 5c 5b 



^ r , L . 6H L 5H JH 5H , [5H\ 

„ , , 5H 5H 5H 5H ( 5H\' 
YH{c) = lb a+ ld C - C Ja-- d lb-{-5b) ' 

To derive (A.8) we have used that for a function H of the form (A. 2) 

Yff(d) = 



-5d' d 



f5H\ 6H, 5H 
{id) + lb b - C lc- = °- 



This follows from the fact that 



for any 



V(0 ) = V(0o) 



(A.9) 



(A.10) 



(Alia) 



e 



e D ]9 ° 



e 



-D 



+ 



e 



D 



e 



-D 



D e gl s . (Allb) 



Relation (A. 11) is easily verified directly, or it may be traced back to the symmetry in 
(3.1b). Writing the arbitrary variation 5H of H in (A. 2) in the alternative ways, 

J trres (([5a - 5b(d + d)~ l c - b(d + d)~ 1 5c + b(d + d) _1 <ta(<9 + d)' 1 ^ jj-(K)) 



(A12) 



leads to the relations 



(A13) 



Using (A.l), (A.8), (A.9), (A. 10) and (A. 13) it is in principle possible to verify (A.6) directly 
for an arbitrary function 7i, but we find it easier to do so for a conveniently chosen complete 
set of functions on A, which is enough by the general properties of the PB. 
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The complete set of functions for which we are now going to verify (A. 6) consists of 
the functions Hk of the form 

H k (L) = Tr(L£d k ), V£;eZ, V£ e gl p , (AAA) 

for which = £d k . In order to simplify the computation (throughout which we shall 
arbitrarily fix £), we make use of the "integrating factor" trick similarly to that in Sect. 4, 
that is we write 

(d + d) = W~ 1 dW, d = W~ 1 W, (A15) 
in analogy to (4.33). We then readily find, 

(£d k K)+ = = (K£d k )+ for k < -1, (A16) 

(^d- 1 K) + = ^=(K^d- 1 ) + for fc=-l, (A17) 

and for k > 0, 

(£d fc if) + = £d fc if + f (M^ -1 )^ -1 ^ = id k K + £D k (b)(d + d) _1 c, 

(A. 18) 

(K£d fc )+ = if£d fc + (-l) fc 6W r - 1 a _1 (W r ce) (fc) = + (-l) fc 6(d + d) _1 Z> fc (cO, 

where we used the covariant derivatives defined similarly to (4.38b), using <i in place of w. 
(Incidentally, w = d in the 6-gauge.) These relations and (A.l) give the l.h.s. of (A. 6) as 
follows: 

X Hfe (K) = for fc<-l, (A.19) 
X H _ X (K) = [£, K] = [£, a] - £' - £&(<9 + d)"^ + b(d + d)- 1 ^, (A.20) 

and for k > 0, 

X Wfe (K) = - = ((-1)*6Z>* (c£) (6)c) - 

(a£P fc (6) + V ({,V k (6))) (9 + d)- x c - b(d + d)- 1 ((-l) k V k+1 (cO - {-l) k V k (c£) a) 
+ 6(<9 + d)- 1 (ciV k (b) - {-l) k V k (c£)b} (d + d)- 1 c. 

(A21) 

We now set Hk '= 'Hk ° V and verify the equality (A. 6). This equality is trivial for 
k < — 1 and it is also clear for k = — 1 by comparing (A.20) and (A. 8-9) for H_\ using that 

<¥-<• 

as follows from (A. 13). To deal with the case k > 0, we note from (A. 13) that 

^=0, ^ = -(-l)^( c £), 6J ^ = -& k (b), forfc>0. (A23) 
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The explicit form of ^J^- will not be needed, only its property (A. 10). Thanks to this and 
(A. 23) we can rewrite the last term of (A. 21), 



(<9 + d) _1 c 



b(d + d)- 1 (civ k {b) - {-i) k v k {ci)b) (d + d)~ l c = b(d + d)- 1 
= b 6 4±(d + d)-ic-b(d + d)-^c. 

od od 

(A24) 

Plugging this back into (A. 21) and using (A. 23) again, we finally obtain 

Sr (r\ ( 5Hk h 5H A^( SH k SH k SH k fSH k \\ x 

/ x (^-25) 

, Mx n-if^ &H k 5H k (6H k \\ _ n 

~ b{d+d) {-w a+ ^d- c - d ^b--{-w) ) for ^°- 

This immediately yields (A. 6) by taking into account (A. 8-9) and that = for k > 0. 
This then completes the proof of Proposition 4.1, whose claim is equivalent to relation (A. 6). 

Appendix B: The formula of the reduced PB in the DS gauge 

The purpose of this appendix is to present the computation leading to formula (4.34) 
of Sect. 4 that describes the reduced current algebra PB on A4 re d — A^ds- We recall that 
the constrained current J G gl pr+s defining a generic point of .Mds has the form 

r 

J = jp,r,s + C+ = ^2 e r,i ® Vr-i+1 + e r , r +i <8> (+ + e r +i,i <8> C- + e r+ljr+1 ®w + C+, (S.la) 



i=l 



where the explicit matrix form of j p , r ,s is given by equation (0.8) with the variables 

Vi G gl p , w G gl s , C+ £ ma *t (p X s) , C- e mat ( s x P) > (B.lb) 

and C + = J2 k= \ e kjk +i ®T is the constant matrix appearing in A Pjr)S = C + + AC_ in (0.6). 
We know already by Theorem 4.2 that the PBs of functions of the components vi, . . . ,v r 
are given by the standard quadratic Gelfand-Dickey PB on the space of operators £, 

r 

£ = L+ = A r d r + J2 u ^ r ~\ u i = Av i A r - i , A = -r" 1 . (B.2) 

i=i 

We need then only to compute the other PB relations. We choose to do this by computing 
the Hamiltonian vector field X# := for H = Q, P, R respectively where 

Q(J) = tr (/C+) , P(J) = r tr (^C-) , R{J) = r tr ( aw ) ( S - 3 ) 
Jo Jo Jo 
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with /, <p, a being matrix valued test functions. It is not hard to see, for instance from the 
theory of reduction by constraints, that Xp takes the form 

£ TT 

X H (J) = [K H , J] — K' H , K H = — + B H , (BA) 

where [-Bp, J] — B' H is a linear combination of the Hamiltonian vector fields associated to 
the ("second class") constraints that define the above special form of J. The polynomial 
nature of the DS gauge ensures (see e.g. [20]) that once J and are given one can uniquely 
solve (B.4) (where the form of Xp(J) must be consistent with that of J in (B.l)) for £?p, 
X// in terms of differential polynomial expressions in the components of J and ^j. 

Let us determine in turn Xq, Xp, Xp. Inspecting the simplest examples leads us to 
search for Kq in the form 



r 



Kq = ^2 e r+l,r+l-i <S> fi 



i=l 



/0 ••• 0\ 

••• 

\fr ■■■ fl 0/ 



with h = f. (5.5) 



Substituting from (B.5) in (BA) we get 

X H = K+, Xq(C+)=0, X (C_) = /t; r -X>(/ P ), X o (t;0 = -C+/i, (B.6) 
and the recursion relation 

fi+i = fviA - V(fi)A, (B.l) 

where T>((3) is given by T>((3) = (3' + w(3 for any (3 G mat(s x p). The solution of this 
recursion relation is found to be 

i-1 

fi = (-l) i ~ 1 V i ~ 1 (f)A i ~ 1 - ^2(-l) i - k V i - 1 -\fv k )A i - k . (B.8) 

k=i 

Substitution from (B.8) in (B.Q) gives Xq explicitly. 

To do the analogous computation for H = P, it is advantageous to change variables by 
transforming to a different gauge section. In fact, there exists a unique gauge transformation 

J~J = gJg- 1 -g , g-\ g=(* ° ), (B.9) 



where A is a block lower triangular matrix with p x p unit matrices along the diagonal, for 
which J takes the form 

r 

J = ^2 e M ® Vi + e r,r+l <8> C+ + e r +i,i ® C- + e r+ i, r+ i <S> w + C+. (B.10) 
i=l 
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The 5j's in (B.10) are unique differential polynomials in the v^s in (B.l), which may be 
determined from the equality 

r r 

A r d r + J- Ul d r ~ l + z + (d + w)-^- =L = A r d r + dr ~ 1 ^ + z+(d + w)- x z- (£.11) 
i=l i=l 

using the notations 



z+ = -AC+, z.=C-, u i = Av i A r ~\ ui = A r+1_i {;j. 



(B.12) 



This equality results from the elimination procedure performed in the respective gauges 
(B.l) and (B.10). In the latter gauge we have a formula for Xp precisely analogous to 
(B.4), 

(B.13) 



X P ( J) = [K P , J] -K' P , Kp = T1 + Bp, 

oJ 



and Kp turns out to have the form 



K P = ^ e i,r+i ® <Pi 



i=l 



/0 ••• i Pl \ 



••• <p r 

Vo ••• o o / 



with <pi = cp. (5.14) 



In fact, substituting from (-B.14) in (5.13) leads to 

X P (w) = -£-<p, X P (C+) = -T)(<p r ) - v r <p, Xp(C-) = 0, Xp(vj) = ipi(- (5.15) 



with the recursion relation 



where V(j3) is given by V(j3) = (3' — (3w for any (3 G mat(p x s). This yields 



(5.16) 



i-l 



= A i - 1 V i ~ 1 { v ) + ^2A i ~ k V i ~ k ~ 1 (v k <p). 



(B.17) 



k=i 



Plugging this back into (B.15) gives Xp explicitly. 

In the case of H — R we find that K R equals |j, i.e., K R = e r+ i yr+ i (g> a. Therefore 

^-r(w) = [a,w] - a', Xp(C+) = -C+«, X fl (C-) = «C-> Xp(^) = 0. (5.18) 



The remaining task is to find a neater form of the above formulae. First we rewrite 
them in terms of the operator £ in (B.2) as follows. 
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Claim Bl: Formula (5.6) with (5.8) is equivalent to 

X Q (£) = (A(+(d + w)- 1 fT£) + , X Q (C+) = 0, Xq(C-) = -V(f r ) + fv r , X Q (w) = f(+. 

(5.19) 

Claim B2: Formula (5.15) with (B.17) is equivalent to 

x P (£) = (e< f (d + w)- 1 c_) + , x P (c + ) = -v( Vr )-v r <f, x P (C-) = o, x P H = -C-v?. 

(5.20) 

Claim B3: Formula (5.18) is equivalent to 

X R (£) = 0, Xp(C + ) = -C+a, X fl (C-) = aC-, X fl («;) = [a, w] - a'. (B.21) 

Claim B3 is obvious from (B.18). Claim Bl may be verified as follows. We use our 
usual trick to write 

8 + w = W~ 1 dW with w = W~ 1 W. (B.22) 

Then we have 

r-l 

-A( + (d + w)- 1 fT£= A( + W- 1 d~ 1 WfA r - 1 d r - A ^ ( + W~ 1 d~ 1 WfTu r - i d i . (5.23) 

Using the identity V{F) = W- l {WF)' and the formula d^F = E^qC -1 )*-^"* -1 we can 
write the contribution to non-negative powers of d in the right hand side of (B.23) in the 
form 

r — 1 r — 1 i—1 

AC+ £(-l) fc Z>*(/)A'- 1 r -*- 1 - AC+ ^(-^(frur^d^- 1 (5.24) 

fc=0 i=l fc=0 

which simplifies to 

r r 

£ AC+/ fc A r - k dT- k = -J2 *Q(u k )d r - k , (5.25) 
fc=i fc=i 

as required by Claim Bl. Similarly, one may check Claim B2 using the identity V (^F^j = 
(FW^yW and the expression of £ = L + in the variables Ui provided by (B.ll) and (B.12). 

Let us now consider the Hamiltonian H = Hi + Q + P + R given by 

p2,7T p2,7T p2,7T /»2"7T 

H(J)= / trres(££) + / tr (C+/) + / tr(C_^) + / tr (wa) , (5.26) 
Jo Jo Jo Jo 

where £ is a p x p matrix PDO of the form £ = $^i=i &c?* -r-1 with arbitrarily chosen 
£i G (7 / p . We wish to present the Hamiltonian vector field X# associated to H by means of 
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the reduced current algebra PB in terms of the variables £, w, and z± in (B.12). To make 
contact with formula (4.34) of Sect. 4, we now substitute 



/ - - 6H A <p - 6H a- 6H £- SH 
5z+ ' 5z- ' 5w ' bi 



(B.27) 



Combining the above claims with Theorem 4.2 of Sect. 4 implies the following formula: 



,6H. 



SH 



SH 



X H (z + )=P (£ 6 -^z + W- 1 



Xh(z-) = -W- 1 ^ I Wz-^-i ) - (-iyv 



W + A r V r 



SH 



5z. 



r—k 



z + (d + w) 
SH 



.i 5# 

5z + 



A:=l 



SH 

Sw ' 



5tf 



Si 



SH 

Sz + 



V(-i) r - fc p r - fc 



fc=i 



5w 



Xh(«j) 



JzZ 



+ 



SH_ 

Sw 



w 



SH_Y 
Sw ) 



(5.28) 

Here the first term in ~X.h(z±) has been found using the skew-symmetry property of the PB. 
Recall that the notations P and P ' are defined as follows. For any PDO %, expand \ so that 
all powers of d are on the right, then Po(x) i s the coefficient of <9°; alternatively writing 
X so that all powers of d are on the left, Pq(x) is the coefficient of d°. This definition 
allows to convert (B.28) into formula (4.34) given in Theorem 4.4, which completes our 
derivation of this formula. The reader may derive (4.35) from (4.34) by computing the Lie 
derivative mentioned in the proof of Theorem 4.4. We only note that, like in the above, it 
is advantageous to make use of the linear functions of J for this purpose. 
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